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I. Experimental Investigation* 


A. von Hippe., F. G. Cuestey, H. S. DeNMaARK,** P. B. ULIN, AND E. S. RITTNER 
Laboratory for Insulation Research, Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received March 4, 1946) 


An experimental study of thallous sulfide photo-conductive cells is presented. It includes an 
investigation of the sensitizing oxidation reaction; cell characteristics as function of voltage, 
temperature, light intensity, wave-length, modulating frequency, and composition; measure- 
ments of thermoeffect and optical absorption; x-ray and electron diffraction studies, and an 
examination of the photo-response to light flashes of very short duration. A theoretical dis- 
cussion will be given as Part II in a subsequent paper. 





INTRODUCTION 


HE photo-sensitivity of oxidized thallous 

sulfide was discovered by T. W. Case in 
1917. Since his initial work! some progress has 
been made in improving the stability and re- 
sponse of this photo-conductive cell** as a result 
of empirically developed*-* and undisclosed? 
processes of manufacture. More recently the 
photo-cell, as further developed by R. J. Cash- 
man,’ has become of importance in this country 
for war communication systems, but its use- 
fulness at present is limited by lack of appreciable 
sensitivity above a wave-length of 1.24 and by 


_ This research was carried out under contract OEMsr- 
1036 with Division 16, Section 16.4, of the National 
Defense Research Committee and completed in June, 
1945, Publication has been delayed because of the war. 

* Deceased. 

. lr. W. Case, Phys. Rev. 15, 289 (1920). 

a R. Sewig, Zeits. f. tech. Physik 11, 269 (1930). 

iW. Leo and C. Mueller, Physik. Zeits. 36, 113 (1935). 
(1928) Majorana and G. Tedesco, Atti accad. Lincei 8, 9 

at Michelssen, Zeits. f. tech. Physik 11, 511 (1930). 

“A. A. Sivkov, J. Tech. Phys. U.S.S.R. 8, 11 (1938). 

R. J. Cashman, N.D.R.C. Report No. 16.4-6 (1943). 


poor frequency response.* A fundamental study 
was therefore undertaken with the objective of 
gaining insight into the photoelectric mechanism, 
as such knowledge should be of considerable 
assistance in understanding and improving the 
cell properties. 


I. PREPARATION OF THALLOUS SULFIDE 


Three batches of thallous sulfide were pre- 
pared* by two different methods. In the first of 
these, thallous sulfide was precipitated by means 
of hydrogen sulfide from an ammoniacal solution 
containing thrice recrystallized thallous sulfate. 
The deposit was washed, dried, and fused in 
vacuum, yielding a product of reasonable purity 
(see analysis of sample No. 1 in Table I). How- 
ever, owing to unavoidable oxidation in the course 
of handling the finely divided precipitate under 
atmospheric conditions, this batch of material 
was used only for exploratory experiments. 


8 W. L. Hole and R. J. Cashman, N.D.R.C. Report No. 
16.4-7 (1943). 

* The preparation of the material was the work of M. C. 
Bloom and J. P. Levy. 
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TABLE I, Spectrographic analyses. 








Impurities, parts per million 





Material Al Cu Fe Ca Ag Na Mg Ni Pb B 
TLS No. 1 (pptn.) 6 4 4 25 — 40 10 5 35 — 
Tl:S No. 2 (synthesis) 0.5 2 0.5 0.5. <.04 04 A <1 > 100 — 
ThLS No. 3 (synthesis) 1 8 <1 1 <.05 <1 <.3 <2 3 <2 








The other batches were prepared by reacting 
stoichiometric quantities of thallium and sulfur 
at elevated temperatures in evacuated quartz 
apparatus. The sulfur was purified by the method 
of Bacon and Fanelli.? In the case of sample 
No. 2, C.P. thallium metal, which had been 
melted and strained in vacuum to remove oxide, 
was employed. A more elaborate procedure was 
followed in purifying the metal in the preparation 
of sample No. 3, consisting of three successive 
fractional crystallizations of thallous sulfate, 
cathodic deposition of thallium from the satu- 
rated sulfate solution acidified with dilute sul- 
furic acid, washing, drying, and melting and 
straining in vacuum. This latter sample was also 
stirred with nitrogen during the reaction to 
assure homogeneity throughout the melt. 

Table I summarizes the spectrographic analyses 
of the different thallous sulfide products. 

Sample No. 2 was employed in the majority of 
the experiments described in this paper. Sample 
No. 3, available only in the later stages of the 
research, served to demonstrate that the photo- 
sensitivity of thallous sulfide originates from its 
oxygen content and not from accidental im- 
purities. 


Il. THERMAL STABILITY OF THALLOUS SULFIDE 


Inasmuch as thallous sulfide is evaporated in 
vacuum as the first step in preparing photo-cells, 
its thermal stability is of great importance. From 
known values of free energy content and heat of 
formation, calculation* was made of the change 
in free energy attending the change in state: 


T1.S (solid)—>2 TI (liquid) 
1 
+ (gas, 10-* mm) at 450°C. 


The calculation assumes that S, is made up of an 


(1942) F. Bacon and R. Fanelli, Ind. Eng. Chem. 34, 1043 
* Calculation performed by J. H. Schulman. 


approximately equal mixture of Ss and Ss mole- 
cules. AF is positive in sign and amounts to about 
20,000 calories, indicating that thallous sulfide 
will not spontaneously decompose into its ele- 
ments under the above mentioned conditions. 

An attempt was made to confirm this con- 
clusion by heating a thallous sulfide layer, 
evaporated in the usual type of cell (see Section 
IV), for several hours in vacuum to progressively 
higher temperatures. After each heating period 
the cell was quenched and measurement made of 
its conductivity, which certainly should be 
affected by decomposition. No appreciable 
change in conductivity resulted until at about 
400°C the experiment came to an end because the 
layer sublimed off the grid electrodes. 


Ill. OXIDATION REACTION 


Thallous sulfide is photo-sensitized by oxida- 
tion. A study of the oxidation reaction at 250°C 
and 350°C by gas absorption, and x-ray and 
electron diffraction methods was therefore under- 
taken to identify the reaction products. 


A. Gas Absorption Method 


It is a prerequisite of accurate gas absorption 
measurements that the system be of small volume 
so as to maximize the change in pressure ac- 
companying the reaction. The essential portion of 
the apparatus employed is sketched in Fig. 1. 
Volume B between the stopcocks C and D was 
calibrated with water before sealing to the 
vacuum system (B~45 cm?). The small volume 
of the Pirani gauge,* P, was estimated from the 
bore of the tubing employed. The reaction vessel 
comprising the major portion of volume A (A ~45 
cm?) was made from Nonex glass tubing of 
30 mm o.d. 

For each reaction vessel, volume A was cali- 


* A description of a more rugged model of this gauge, 
which is capable of measuring pressures up to 15 mm, 
will be published elsewhere. 
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brated by filling A and B to a predetermined 
pressure, evacuating B, reconnecting A and B, 
and establishing the resulting pressure. A weighed 
amount of thallous sulfide (14.5+1 mg in the 
usual case, known to +0.1 mg) was then intro- 
duced into the reaction vessel, the system 
evacuated, the vessel outgassed at 400°C for an 
hour, and the entire sample evaporated by 
heating with a Bunsen burner (film thickness 
about 5000A). To initiate the reaction, the vessel 
was heated by an oven to the desired tempera- 
ture, tap D closed, oxygen at known temperature 
introduced into volume B to the desired starting 
pressure, and tap D opened. The progress of the 
reaction was measured by observing the pressure 
as a function of time at the constant temperature 
for a period up to twenty hours. The reaction 
vessel was then quenched by an air blast and the 
pressure re-measured at room temperature with 
the whole system in temperature equilibrium. 
This completes the data required to calculate the 
composition of the oxidized layer. 

Visual observation of thallous sulfide layers 
oxidized to varying degrees reveals that the 
photo-cell material consists of two phases: a 
grayish-black (photo-sensitive) substance and a 
bright yellow (insulating) material.. With pro- 
gressing oxidation the quantity of the yellow 
phase increases; and at a composition approxi- 
mated by 1 mole O2 per mole T1.S, the black 
phase completely disappears. 

The existence of these phases, and of only these 
two phases, was confirmed by the x-ray and 
electron diffraction studies reported below. 


B. X-Ray Diffraction Study 


X-ray powder photographs were taken of a 
number of photo-sensitive layers carefully re- 
moved from the cells in a nitrogen atmosphere 
and sealed in small glass capillaries. A phase 
found to be present in all sensitive cells gave rise 
to a diffraction pattern* correlating satisfactorily 
with the structure of T1,S proposed by Ketelaar 
and Gorter.!° Usually a second (yellow) phase 
was observed in addition. The diagram in Fig. 2 
shows the positions and relative intensities of 


tintin 


* This same pattern was obtained with our fused samples 
of thallous sulfide and with unoxidized evaporated layers. 

© J. A. A. Ketelaar and E. W. Gorter, Zeits. f. Krist. 
101, 367 (1939), 


diffraction lines from patterns obtained with pure 
T1.S, with a typical photo-cell sample and with 
the yellow material from a ‘‘completely”’ oxidized 
layer. 

Inasmuch as the photo-sensitive phase pos- 
sesses the same structure as thallous sulfide, it is 
probably safe to conclude that the former repre- 
sents a solid solution of oxygen in T1.S. The 
change in lattice constants with composition 
usually observed in solid solutions, if present, 
could not be detected in the present case because 
the small crystallites in the evaporated layers 
produced diffuse lines in the high @ region and 
prevented accurate measurements. 

The interplanar spacings calculated from the 
diffraction photograph of the yellow phase indi- 
cate a face-centered cubic lattice with a cell edge 
of 10.77A. Since this phase has a structure 
differing from any of the known structures of 
thallium compounds and since its composition is 
about 1 mole O2/mole T1.S, it may be identified 
as the compound T1,SO>."! This substance proved 
quite stable; it cannot be decomposed at a 
reasonably rapid rate by heating in vacuum 
unless a temperature high enough to sublime the 
material is employed. The gray condensate re- 
sulting consists of a mixture of the Tl,S and 
T1,SO2 phases. 

In passing, it may be noted that no changes in 
structure were observed on heating thallous 
sulfide up to 385°C. From the displacement of the 


REACTION VESSEL 


cre 


Fic. 1. Reaction vessel, system A. 














1! Reference to the existence of this compound has been 
made by I. I. Iskoldsky, Mineralnoje Syrye No. 4, 404 
(1931). 
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Fic, 2. X-ray diffraction patterns. 


(0006) reflection in the x-ray diffraction patterns 
it was possible to determine the linear coefficient 
of thermal expansion in the direction of the ‘‘c”’ 
axis. The value, a=34X10-* (20-385°C), does 
not differ appreciably from that of thallium metal, 
a = 30.3 10-8. 


C. Electron Diffraction Study 


For the electron diffraction investigation a cell 
was produced in the diffraction apparatus itself 
for a study at grazing incidence. On account of 
the relatively low conductivity of the thallous 
sulfide, the diffraction specimen charged up if not 
prepared on a conducting surface. The sample 
prepared on a glass slide consisted of two parts: a 
section on aquadag from which the diffraction 
studies were made and immediately adjacent to 
it a photo-cell section with its standard gridwork 
(see Section IV). Thallous sulfide was vacuum 
evaporated from a small crucible onto the surface 
of this slide and also onto Formvar Films of 
thickness about 100A for transmission studies. 
Examination of the diffraction patterns obtained 
(see*Fig. 3, a+-transmission patterns, c-reflec- 
tion pattern) show the crystals of thallium 
sulfide to be arranged with their “‘c’’ axis pre- 
ferably oriented in a direction perpendicular to 
the receiving surface. The slide was then heated 
for twenty minutes at 250°C at an oxygen pres- 
sure of 2 mm, resulting in the development of an 
appreciable photo-sensitivity. The electron dif- 
fraction pattern from the sensitive surface is 
shown in Fig. 4a; the interplanar spacings check 
those from the unoxidized T1.S pattern. In addi- 
tion, it may be noted that the original orientation 


VON HIPPEL, CHESLEY, DENMARK, 


ULIN, AND RITTNER 


persists and that the crystallites have grown in 
size. However, the pattern is considerably more 
diffuse than that from the unoxidized surface and 
indicates that the structure has probably been 
disrupted by solid solution of some oxygen. More 
complete oxidation produced the pattern shown 
in Fig. 4b; the interplanar spacings check the 
x-ray data from the yellow material. Hence, the 
electron diffraction data lead to the same con- 
clusions as those suggested by the x-ray diffrac- 
tion results. 


IV. PREPARATION OF CELLS 


Cells prepared by a single oxidation treatment 
at 250° and 350°C were investigated with the 
purpose of obtaining data on dark resistance and 
photo-sensitivity as a function of oxygen content. 
In addition, the reproducibility of the charac- 
teristics of photo-cells prepared under closely 
identical conditions was studied. 


A. Methods 


Two vacuum systems of different design were 
employed in the preparation of these cells. The 
first one, system A, has already been described in 
Section III A and was used with the following 
changes. The reaction vessel was furnished with 
sealed-in tungsten electrodes and an aquadag 
gridwork as shown in Fig. 5.* The gridwork 
covered an effective area of 1X1" with 16 lines 
per inch. Pretreatment of the cell blank consisted 
of a heat treatment in air for several minutes at 
400°C followed by a heating in a vacuum for one 
hour at 500°C. The introduction of sample, out- 
gassing, evaporation, and oxidation were carried 
out by a procedure similar to that outlined in 
Section II] A. After quenching, the resistance and 
photo-sensitivity of the cell were determined ; the 
oxygen was pumped out, which did not affect the 
sensitivity ; and the cell was sealed off. Since its 
properties changed rapidly with time for a few 
days, the resistance and sensitivity increasing in 
general, final measurements were postponed until 
the characteristics had stabilized. 

The essential features of the second system, 
system B, which was designed to operate with a 
large volume (2 liters) of oxygen at pressures of 
the order of those employed in the commercial 


* After a design of R. J. Cashman. 
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(c) 


Fic. 3. Electron diffraction patterns of thallous sulfide: 
(a) electron beam | to T1.S film a 
(b) Tl2S film inclined 45° to electron beam /t™@™smission 
(c) electron beam at grazing incidence to T1,S surface. 


Process,” are shown in Fig. 6. The preparation of exceptions that the volume of the cell blank, 


cells was carried through as before with the minor ca. 40 cc, was not accurately determined, that 
A 


, ras i i eating, and tha 
Dae pHewlett. J. J. FitzPatrick, H. T. Wrobel, OXYSen was introduced prior to heating, and that 
N.D.R.C. Report 16.4-39 (1945). the reaction was initiated by covering the cell 
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(b) 


Fic. 4. Electron diffraction patterns of oxidized thallous sulfide layers: (a) Tl:S phase, (b) TleSOz2 phase. 


with the hot oven. The oxidation was allowed to 
proceed for about twenty hours; at the end of this 
period the cell was quenched, sealed, and meas- 
ured as described above. The residual pressure 
was also determined to permit an estimation of 
oxygen content.* 


B. Resistance and Photo-Sensitivity as Function 
of Oxidation 

Data on dark resistance, thermal activation 

energy, U, and photo-sensitivity of cells, pre- 


30 MM NONEX GLASS TUBING 





040" W WIRE 





L | 6M NONEX 
’ GLASS TUBING 
AQUADAG GRIDWORK 


bets +434 





Fic. 5. Cell blank. 


* We are indebted to R. B. Speed for the preparation of 
cells in system B. 


pared from sample No. 2 in systems A and B by 
oxidation at 250°C, as a function of the composi- 
tion of the oxidized layer are tabulated in 
Table II and similar data at 350°C for system B 
are shown in Table III.** 

For the cells prepared at 250°C there is no 
discernable relation between thermal activation 
energy (see Eq. (4)) and oxygen content; how- 
ever, the dark resistance increases with oxidation 
until quite sensitive cells are obtained. Beyond 
this point correlation between resistance and 
composition is obviated by the large fluctuations 
in the value of the dark resistance of cells pre- 
pared by the same procedure. 

Table II shows, in addition, that the minimum 
amount of oxygen required for producing a 
sensitive cell varies with the range of pressures 
employed. Although under the conditions of 
system A about 0.3 mole O2/mole T1,S is re- 
quired, in system B about 0.1 mole O2/mole Tl:S 
suffices. This apparently results from the simul- 
taneous or successive occurrence of two oxidation 
products, the relative amounts of which vary in 
accordance with operating conditions. The ab- 
sorption of still more oxygen does not improve 


** See Section V for details on measurements, definitions 
of terms, and references. 
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the photo-conductivity further but converts more 
of the thallous sulfide phase into Tl,SO2. Highly 
sensitive cells can therefore be made at 250°C 
even though the operating pressure varies over 
wide limits. The variations in sensitivity of cells 
prepared by the same process are appreciable but 
not as extreme as the fluctuations in the resistance. 

In the case of the cells prepared at 350°C, the 
compositions are higher in oxygen than the 
corresponding ones attained at 250°C owing to 
the increased rate of reaction at the higher 
temperature. The electrical characteristics do not 
differ fundamentally, however, in the two sets of 
cells. 


V. CHARACTERISTICS OF THALLOUS SULFIDE 
PHOTO-CONDUCTIVE CELLS 


A determination of whether the photo-sensi- 
tivity of the cell represents a primary or second- 
ary photo-effect is of paramount importance. 
The sluggishness of the cell," its strong 
temperature dependence’ and poor frequency 
response? indicate a secondary type. In addition, 
Runge and Sewig'® have concluded from their 


rl) 








= 
Sy 
Fic. 6. Reaction vessel, system B. 


measurements that light creates disturbance 
centers which decay exponentially with time. An 
unequivocal decision can be obtained by a 
determination of the quantum yield. While a 
primary photo-effect makes use only of the 
charge carriers liberated by the absorption of 
light quanta and is normally limited to a maxi- 
mum yield of one electron per quantum absorbed, 
a secondary photo-effect changes the conductivit y 
of the material and quantum yields far exceeding 
unity may result. 

The characteristics of representative Cashman 
cells* were investigated as a function of voltage, 
temperature, light intensity, and wave-length. 


TABLE II. Resistance and photo-sensitivity vs. oxidation, 250°C, 








Composition, Initial 
mole O:2/ pressure, 
mole T1l:S mm 


Dark resistance, electron 


White light sensitivity 
Response to 90- 
cycle modulated 

activation light, signal: 

energy, AI /Io at noise with infra- 

.023 hololumens red transmitting 

megohms volts per in.? filters db 


Thermal 





Qb 

0.119 
0.251 
0.339 
0.578 
0.760 
14 


.06 
13 
.16 


39 
49 


.0075-.022 59 — 
O11 48 d 0 

6.1 62 . 27 

28 A7 66 
225 61 40 
45 46 57 


oe _— — 
.022-.11 54° — 

0.55 52 4 24 
187 79 62 
1.7-204 77° 19-198 52-62 


2.5-55 . 18-44 58-67 
37 41 35 61 








* Without filter S/N ratios are about 12 db higher. 
> After initial evaporation. 


° Owing to the high conductivity of the unoxidized material, the resistance of the gridwork is not negligible compared to that of a layer of standard 


thickness; hence, a much thinner layer was employed. 
4 Sample completely yellow. 
® Cell No. 64. 
f Cell No. 54. 


®W.W. Coblentz, Bur. Stand. J. Research 16, 254 (1920). 


*T. Kaplan, J. Opt. Soc. Am. 14, 253 (1927). 
[T. Runge and R. Sewig, Zeits. f. Physik 62, 726 (1930). 


* The gridwork covered an effective area of ?’’X ?” with a spacing of 1/30’’ between adjacent lines. 
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TABLE III. Resistance and photo-sensitivity vs. oxidation, 350°C. 








Composition, 
mole O2/ 


Sample Cell No. mole Tl2S mm 


Initial 
pressure, 


White light sensitivity 
Response to 90- 
’ cycle modulated 
light, signal: 
noise with infra- 
red transmitting 
filter, db 


AI /Io at .023 
hololumens /in.? 


Dark resistance, 
megohms 





0.21 

0.30 
0.36 
0.38 
0.16 
0.69 


0.10 
0.18 
0.18 
0.16 


138 
168 
166 
136 
116 


114 0.53 


0.066 


72.5 
10.5 
140 
1.05 
0.15 
840 


48 
54 


bdo Ww bdo 
iD Go im 
oo oo Un OO SION 








* Material employed in the commercial process. 


A. Apparatus 


Monochromatic** radiation was obtained by 
the use of a Gaertner instrument (L234-HT) with 
quartz optics, the output of which was focused 
by a mirror system alternately on the test cell or 
on a thermocouple. 

For determining its temperature dependence, 
the cell was mounted in an electrically heated 
housing or in a double-walled container thermo- 
stated by circulating liquids. 

Three tungsten lamps, covering a wide range 
of light intensities and calibrated in terms of 
microwatts as well as microhololumens*** were 
employed in ‘‘white light’’ measurements. Tests 
with square-wave modulated light of one 
microhololumen peak (without filter) were car- 











os 





Za 























225 45 675 90 
VOLTS 


Fic. 7. Dark current, Jo and total current, J, obtained 
with illumination as function of applied voltage (A =0.8y; 
F=0.5 microwatts). 


** AN~I150A. 

*** One hololumen is that radiant flux of all wave- 
lengths which has the spectral distribution characteristic 
of a tungsten lamp at color temperature 2848°K, and which, 
evaluated as visible light, equals one lumen. (See N.D.R.C. 
Report No. 16.4-10 (1943) by G. A. Van Lear, Jr.) 


ried out on the standard instrument developed at 
the University of Michigan.'®!” 


B. Effect of Applied Voltage 


For a thallous sulfide cell the dark current, Jo, 
and the current, J, obtained under illumination 
are proportional to the voltage applied, V, at 
constant temperature (Fig. 7); that is 


Ih= VGo= V/Ro, (1) 
I=I,)+AlI=VG 
V V 


= V(Go+ AG) =—= —, (2) 
R,—AR 


where R> is the dark resistance, R the resistance 
under illumination, and Gp and G the correspond- 
ing conductances. 

It follows from Eqs. (1) and (2) that the cell 


24; 
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. 3 ee 0 
WAVELENGTH, MIGRONS 

Fic. 8. Effect of spot size on photo-current for con- 
stant light input. (Ro=5.4 meg.; V=22.5 v; F=1 micro- 
watt.) 


( if L. Hole, N.D.R.C. Report No. 16.4—Special No. 9 
1944), 

17L. N. Holland, N.D.R.C. Report No. 16.4—Special 
No. 6 (1944). 
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Fic. 9. Temperature dependence of the dark conductance. 


sensitivity can be characterized independently of 
V by the ratios 


AI/Ip=AG/Go=AR/R. (3) 


In addition, presentation of data in this form 
permits comparisons of sensitivity between cells 
of different area and grid spacing. 

In order to utilize the above equations, it is 
essential that the photo-sensitive area of the cell 
be uniformly and completely illuminated. These 
conditions have only been approximated in 
measurements requiring the use of the mono- 
chromator by employing a spot size of diameter 
equal to the length of the side of the square 
sensitive surface. The effect of spot size on cell 
response for constant light input is illustrated by 
Fig. 8. 


C. Dark Current and White Light Response 


The temperature dependence of the dark con- 
ductance of oxidized thallous sulfide is typical of 


that of a semiconductor, Go decreasing expo- 
nentially with the reciprocal of the absolute 
temperature over a wide range (Fig. 9, Eq. (4)). 


Go=aeUl*T, (4) 


a=constant, U=activation energy. 

Owing to the high temperature coefficient of 
the dark current, the accuracy of light response 
measurements is limited at very low and very 
high light intensities; but, as is evident from 
Fig. 10, only at very low levels of illumination 
does the photo-current increase proportionally to 
the light input. At higher light intensities the 
response deviates considerably from linearity. 
This behavior explains the previously noted 
observation (Fig. 8) that the output decreases 
when the light is concentrated on a smaller spot. 


D. Spectral Response and Quantum Yield 


The spectral response characteristic of the 
thallous sulfide cell (Fig. 11) varies with the 
illumination level because of the non-linearity of 
the cell. In addition to the photo-current, AJ, the 
quantum yield, 7, is plotted in Fig. 11, assuming 
that every incident light quantum is absorbed 
and contributes to the photo-effect : 


photo-electrons (1.24)(AJ) 





incident quanta 7 (A) (F) 


(AJ = photo-current in » amp., F=light input in uz 
watts, \=wave-length in microns). Because of 
this assumption, the calculated quantum yield 
will be lower than the true yield; nevertheless, 
the former surpasses unity by a large margin. 


ILLUMINATION, MICROWATTS/INAICOLOR TEMP =2850°K) 


Fic. 10. White light response as a function of illumination. 
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Fic. 11. Photo-current and quantum yield as function 
of wave-length for two light intensities. (Ro=2.5 meg.; 
V =22.5 v.) 


This result establishes beyond question that the 
photo-currents observed in thallous sulfide cells 
are the result of a secondary photo-effect. 

Plots of light response vs. intensity for different 
wave-lengths (Fig. 12) are similar in shape and, in 
fact, can be combined into a single curve by 
multiplying the yw-watt scale of each plot by a 
characteristic factor, a. For example, in Fig. 13 
all curves have been converted into the 0.9u 
characteristic by employing the a factors tabu- 
lated in the figure. It appears that a plot of the 
. different a factors against wave-length represents 
a relative spectral sensitivity curve for constant 
light current output. 


NuAp 
a= (6) 


No.9\0.9 


This curve also corresponds to the relative 
spectral sensitivity curve that would be obtained 
at very low illumination levels within the 
linearity range. 

Similarly, all plots of light response vs. intensity 
for different temperatures and for a given wave- 
length (Fig. 14) can be combined into a single one 
by multiplying the y-watt scale by factors b 
(Fig. 15). Again a plot of 6 against temperature 
indicates how the response varies with tempera- 
ture for a constant light current output. 


VI. THERMOELECTRIC STUDIES ON 
THALLOUS SULFIDE 


The type of electronic conduction prevailing in 
semiconductors can be determined from the sign 


DENMARK, 


ULIN, AND RITTNER 

of the thermal e.m.f. Since the mobile charge 
carriers stream from the hot to the cold junction, 
excess conductivity charges the cold junction 
negative, while the migration of holes produces 


the opposite polarity. 


A. Apparatus and Procedure 


Measurements* of the thermovoltage coeffi- 
cient were made on a number of thin films of 
thallous sulfide in various stages of oxidation. 
Figure 16 shows the special type of cell used in 
these measurements. The two tungsten electrodes 
were ground flat and sealed into the Pyrex bulb 
through glass arms, which could be filled during 
the experiments with liquids of the desired tem- 
perature difference (~60°C). The faces of the 
electrodes were covered with aquadag, which was 
subjected to the standard outgassing procedure 
before thallous sulfide was evaporated onto the 
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12. Photo-response vs. light intensity and wave- 
length. (Ro=2.5 meg.; V=22.5 v.) 


* The thermoelectric study was begun and the majority 
of these measurements were made by G. Oster, the re- 
mainder by R. B. Speed. 
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Fic. 13. Conversion of wave-length curves into 
the 0.9 characteristic. 


upper half of the glass bulb. Thermovoltages were 
measured by means of a feedback amplifier of 
high input impedance.'’ 


B. Results 


Thallous sulfide sample No. 1, which was 
known to be slightly pre-oxidized (Section 1), 
exhibited defect conductivity after evaporation 
and the sign of the thermal e.m.f. did not change 
on further oxidation. However, sample No. 2, 
which was completely free of oxygen, manifested 
excess conductivity after evaporation. With in- 
creasing oxidation of this sample the conductivity 
fell, the photo-sensitivity rose, and the thermo- 
voltage went through zero and reversed sign as is 
shown in Fig. 17. Defect conductivity is repre- 
sented in the figure by a positive sign of the 
thermal-e.m.f. coefficient. 

A minimum in the conductivity characteristic 
did not accompany the change in sign of the 
thermo-e.m.f., in contradistinction to a case like 
PbS!* where the material is an insulator at the 
stoichiometric composition and a defect or excess 
conductor in the presence of a defect or excess of 
the metallic constituent. Hence, if sample No. 2 
does not contain an excess of either of its 
constituent elements, the data seem to indicate 
that thallous sulfide is both an excess and defect 
conductor and that oxidation reduces the origi- 
nally dominating excess conductivity to a higher 
degree than the hole conductivity, thus producing 
a reversal of sign. 


'8S. Roberts, Rev. Sci. Inst. 10, 181 (1939). 
' H, Hintenberger, Zeits. f. Physik 119, 1 (1942). 
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Fic. 14. Effect of temperature on photo-response. 
(A=0.8 uw; V=22.5 v.) 


Illumination of an oxidized layer resulted in a 
slight increase in the magnitude of the positive 
thermovoltage, indicating that light increases the 
defect conductivity of the material. 

Thermal e.m.f. measurements have previously 
been made by Hochberg and Sominski?® on 
thallium sulfide samples subjected to a variety of 
heat treatments in air. We agree with these 
authors that high resistant Tl.S shows defect and 
low resistant material shows excess conductivity ; 
but their conclusion that pure Tl.S is a hole 
conductor, changing to an excess conductor on 
oxidation, is contrary to our findings. 


VIL. TRANSMISSION AND REFLECTANCE MEAS- 
UREMENTS ON THALLOUS SULFIDE LAYERS 


Transmission measurements between 3000 and 
8500A* were made on several evaporated films of 
thallous sulfide prepared in cells of normal con- 
struction but lacking gridworks. Typical data 
obtained on very thin films are shown in Fig. 18. 
It is apparent that the absorption increases at 
shorter wave-lengths. 

Measurements of transmission and reflectance 
between 4000 and 7000A** were made using a 
special absorption cell, fabricated from ‘Preci- 
sion Bore” tubing of square cross section and 
containing a Tl.S layer of standard thickness 
(5000A). The transmission proved to be less than 


20B. M. Hochberg and M. S. Sominski, Physik. Zeits. 
Sowjetunion 13, 212 (1938). 

* A Hilger Spekker photometer was employed in these 
measurements. 

** A Hardy color analyzer was used in this work, 
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Fic. 15. Conversion of temperature curves into 
the 27°C characteristic. 


0.1 percent and the reflectance 25+1 percent 
over this spectral range. 

From some additional data obtained between 
8000 and 12,000A*** on layers of standard thick- 
ness, it appears that there is a peak in the 
transmission curve at about 1.0u and that about 
10 percent of the incident light is transmitted at 
this wave-length. 

Over the entire spectral range investigated the 
presence of oxygen does not seem to alter ma- 
terially the shape of the transmission curve. 


VIII. RESPONSE TO FLASHES OF 
SHORT DURATION 


In order to determine the build-up and decay 
characteristics of the photo-current, an investi- 
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Fic. 16. Cell for thermo-e.m.f. determination. 


*** The Gaertner monochromator was employed for 
this purpose. 


INCREASING OXIDATION —~—=- 
TIME OF HEATING AT 250°C, 
MINUTES. 
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Fic. 17. Relation. between dark conductivity, photo- 
response, and thermo-e.m.f. coefficient with increasing 
oxidation. 


gation was made of cell response to light flashes 
of extremely short duration. 


A. Equipment and Procedure 


The incandescent lamp source mounted in 
front of the monochromator slit was replaced by a 
discharge tube of the Edgerton type which could 
be flashed for approximately 3 x 10~* second by a 
condenser discharge. 

The photo-cell was connected in series with a 
battery and load resistance in the usual manner. 
The load resistance was coupled to a cathode 
follower unit with low input capacitance, ap- 
proximately 5 yuf, to eliminate frequency dis- 
tortion, and the cathode follower output was 
amplified by a single-stage a.c. amplifier.* The 
output of this amplifier was further amplified by 
the Y axis amplifier of a DuMont Type 208 
oscilloscope, and the trace, which appeared on 
the screen of a 5LP1 DuMont cathode-ray tube, 
was photographically recorded. In order to ex- 
pand the time scale, a multiple sweep arrange- 
ment was used which was a modification of the 
circuit described by Richter.”! 


* The electronic circuits were constructed by M. Haugen. 
21 W. Richter, Electronics 17, No. 9, 128 (1944). 
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Fic. 18. Transmission of thallous sulfide layers. 


For our analysis it was desired to determine the 
change in conductance as a function of time. In 
order to simplify the interpretation of the various 
photographed patterns, it was necessary to select 
load resistance values, Rr, so as to obtain 
proportionality between output voltage (across 
R,), ME, and change in cell conductance, AG. 
Since it can be shown that 


R,AG(Ro—AR) Io 


AE=R,Al= ’ 
Go(Ro—AR+Rz) 


(7) 





it is apparent that in order to obtain propor- 
tionality between AE and AG it is necessary that 
Rr«R,—AR. This condition was fulfilled in the 
majority of our experiments by employing a load 
resistance of 10,000 ohms. At extremely low light 
intensities, however, it was possible to use a load 
resistance as high as 200,000 ohms. 

The flash response of the thallous sulfide 
photo-cell was studied as a function of voltage, 
wave-length, light intensity, and temperature. In 
addition, cells of different resistance and oxida- 
tion states were compared under closely similar 
conditions. 

The intensity of the flashes was controlled by 
altering the slit widths of the monochromator. 
The intensity of a single flash could not be 
accurately duplicated, however, since the dis- 
charge was not completely reproducible in in- 
tensity and location. Nevertheless, a measure of 
the intensity could be had by determining the 
peak value of AG/G» from an examination of the 
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Fic. 19. Oscillograph trace of response to flash. 


photographed pattern and from the known dark 
conductance. The relative intensities of the 
flashes were varied over a range of the order of 
3000: 1, the lowest intensity being limited only by 
cell and circuit noise. 

For the temperature measurements the elec- 
trically heated oven was mounted inside the 
same shielded container as that enclosing the 
cathode follower and amplifier circuits. The small 
capacitance of the cell, about 10 uyf, increased 
only by about 50 percent in the oven mount. 


——HIGH INTENSITY 
——-— LOW INTENSITY 


4 
TIME, MILLISECONDS 


Fic. 20. Decay characteristics for two light intensities. 
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Fic. 21. Decay characteristic as function of temperature. 


B. Results 


All cells thus far investigated exhibit a response 
characteristic typified by Fig. 19. After the 
occurrence of the light flash, the photo-current 
builds up in about 5 X 10-5 second to its maximum 
value ; its decay with time is exponential save for 
the very beginning of the curve. This can be 
readily seen in a plot of the photo-response on a 
logarithmic scale against time (Fig. 20), where 
the slope is steepest initially and rapidly de- 
creases to a constant value, x. 

The voltage applied to the cell and a change in 
wave-length of the monochromatic light flash 
between 0.6 and 1.1 have no effect upon the 
shape of the decay characteristic. 


Light intensity and temperature, on the other 


TABLE IV. Peak flash response vs. temperature. 














Peak AG/Go ie 
1.45 25 
1.17 35 
1.05 48 
0.67 63 
0.24 76 
0.11 85 


RITTNER 


ULIN, AND 


TIME CONSTANT, (SEC) 103 





37 39 


(7x) 104 


Fic. 22. Temperature dependence of the time constant. 


hand, have a pronounced effect. With increasing 
intensity the initial slope of the log plot becomes 
progressively steeper while the slope of the linear 
portion of the curve remains the same (Fig. 20) 
within the limits of error of about +10 percent. 
With increasing temperature, although the time 
of rise remains unchanged, the cell decay becomes 
progressively more rapid and the peak in the 
response curve falls (Fig. 21). The variation of 
peak response with temperature is a_ linear 
function as the data of Table IV show. The time 
constant, 1/x, of the exponential portion of the 
decay curve increases exponentially with the 
reciprocal of the absolute temperature (Fig. 22). 
In the one instance where reliable measurements 
are available, the value of the activation energy 
(0.65 electron volt) obtained from the slope of this 
curve is slightly larger than the activation energy, 
U (0.54 electron volt), determining the dark 
current dependence on temperature. 

On comparing sensitive cells of different. re- 
sistances and oxidation states (Table V), it 
appears that the time constant, 1/x, increases in a 
very general way with decreasing cell resistance, 
although no correlation between time constant 
and composition can be seen. No relationship has 
been found between the magnitude of the peak 
in the response curve and 1/x, Ro, or composition; 
the data are not included in the table as all 
measurements were not made at the same in- 
tensity level. 

On comparing insufficiently oxidized cells at 
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the same light intensity and applied voltage 
(Table VI—cell No. 4 of high sensitivity included 
as a typical end member of the series), it can be 
seen that the peak yield increases rapidly with 
oxygen content while the time constant remains 


TABLE V. Time constant vs. dark resistance and _ , 
composition of sensitive cells. 











Composition, 


(1/x) 108, sec. Ro, megohms* mole O2/mole T12S 


SULFIDE 


CELLS 369 


TABLE VI. Peak yield, time constant, dark resistance vs. 
composition of insensitive cells. 








Composition, 


(1/x) 10°, Ro, mole O2/ 


Cell No. 


840 
450 
322 
225 
140 
72.5 


_ 
SB te inn 
SoAw } 


0.69 
0.578 
0.307 
0.16 
0.36 
0.21 


Cell No. 


Peak AG/Go 


sec. 


megohms 


mole Tl:S 





59A 
SA 
4A 


1.5 X10-5 
3.23 X10-3 
3.79 X107! 


1.7 
1.2 
2.6 


0.018 — 
0.132 
51 


0 
0.119 
0.25 





Cell No. 


18A 





Frequency, 


1/x, sec.* 


cycles /sec. 


50B 3.6 X10~2 Signal,** db 
3.3 X10~* 


Signal, db 








30 
104 
58 


90 180 


98 93 
5858 


360 


86 
56 








51 0.25 
0.339 

0.39 

0.16 

0.39 

0.30 

0.16 

0.39 

0.39 

OS 0.38 
116B 37 AS 0.16 


UMP ONNSOOCO- 


ReRWUAD 


* The differences in values of Ro from those reported in Table II are 
the result of long aging. 


relatively unchanged. However, the time con- 
stants of cells 59A and 5A are considerably 
smaller than Table V would lead one to expect 
for cells of such low resistance. 


* From flash response data. 
** Without filter; reference level 0.1 volt =0 db. 


IX. FREQUENCY RESPONSE 


The relative response to modulated light as 
function of frequency varies widely for different 
photo-cells as indicated by the data of Table VII. 
Cells possessing very short time constants, which 
are not common, exhibit good relative response 
with frequency over the range investigated. The 
quantitative relation between frequency response 
and time constant will be given in connection 
with a general theoretical discussion to be 
published in a subsequent paper. 
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Thallous Sulfide Photo-Conductive Cells 


II. Theoretical Discussion* 


A. VON HIpPeL AND E. S. RITTNER 
Laboratory for Insulation Research, Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received March 4, 1946) 


A theory is formulated which is in substantial agreement with all facts reported in a previous 
paper on thallous sulfide photo-conductive cells. The underlying physical picture is that light is 
absorbed in the T1.S proper and in the presence of oxygen, positive “‘holes’’ and negative oxygen 
ions are formed. The “holes” and ions recombine at a rate proportional to the product of their 
concentrations. During their lifetime, the negative ions act as space charge compensators for the 
defect conduction and allow correspondingly higher currents to flow. This space charge compen- 
sating action, which is similar to the role played by positive ions in gas discharges, seems to be 
one essential mechanism of producing secondary photo-effects. 











INTRODUCTION 


N a previous paper! an account was given of an 

experimental study of thallous sulfide photo- 
conductive cells. A satisfactory theory should 
account for the following facts established by 
this study: 


(1) the dark current and photo-current are proportional 
to the applied voltage; 

(2) the dark conductance decreases exponentially with 
the reciprocal of the absolute temperature; 

(3) the photo-current, which represents a secondary 
photo-effect, is a non-linear function of the illumination and 
decreases rapidly with increasing temperature; 

(4) the introduction of solute oxygen atoms into the 
TIS phase reduces the conductivity greatly, changes the 
material from an excess to a defect conductor, increases the 
photo-sensitivity many fold, but does not alter materially 
the transmission of the layer; and 

(5) the photo-current produced by a light flash of very 
short duration builds up to a maximum in about 5X 107° 
second and decays exponentially with time except initially 
when the rate of decay is highest and increases with light 
intensity; the time constant of the exponential portion of 
the decay characteristic increases exponentially with the 
reciprocal of the absolute temperature. 


I. THE CONDUCTION MECHANISM IN 
THALLOUS SULFIDE 


Crystalline thallous sulfide, according to Kete- 
laar and Gorter,? consists of TI-S-Tl layers 





* This research was carried out under contract OEMsr- 
1036 with Division 16, Section 16.4 of the National De- 
fense Research Committee and completed in June, 1945. 
Publication has been delayed because of the war. 

1A. von Hippel, F. G. Chesley, H. S. Denmark, P. B. 
Ulin, and E. S. Rittner, J. Chem. Phys. 14, 355 (1946). 

2]. A. A. Ketelaar and E. W. Gorter, Zeits. f. Krist. 
101, 367 (1939). 


arranged similarly to those in the cadmium 
iodide (C6) type of structure. The binding 
between thallium and sulfur atoms cannot be 
purely ionic as this would require the close 
juxtaposition of layers of positively charged 
thallium ions. Moreover, the thallium and sulfur 
atoms can best be accommodated in the structure 
by using atomic radii.2 Similarly, the binding 
cannot be predominantly metallic because the 
electronic conductivity of thallous sulfide is too 
small (10-°-10-* ohm cm at 25°C) for the 
metallic bond to contribute appreciably to the 
cohesion of the crystal. It appears likely, there- 
fore, that the bonding between the sulfur and 
thallium atoms is covalent (possibly with partial 
ionic character) and that the attraction between 
adjacent thallium layers is the result of van der 
Waals forces similar to the cadmium iodide. case. 
The electrical conductivity is probably restricted 
to the TI-TI planes, in which some of the thallium 
atoms lie nearly as close together as in thallium 
metal. 

If the electrons or holes needed in the con- 
duction process were not originally present in 
the compound but were to migrate in from the 
electrodes, space charge formation would restrict 
the current density to very small values. This 
follows from Poisson’s equation : 

dE Arnpe . 
—= (1) 


dx —_ 





where e=electron charge, E=field strength, 
p=density of charge carriers (is a function only 
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of x), and ¢€,=optical dielectric constant. By 
replacing p by the current density, j, 


j=pebE=cE, (2) 


where o=conductivity and 6=mobility; and 
integrating, the following expression results: 


8a jx } 
E= ( +24*) : (3) 


Ex 


Maximum space charge and current density will 
be attained when the field gradient, Eo, at the 
electrode from which the charge carriers emi- 
grate (x=0) is reduced to zero. By introducing 
the voltage, V, applied across the sample of 


length, d, 
d 


v={ Edx; (4) 
0 


and integrating, the maximum current density 
possible under space charge conditions is found: 


V dex 


Jmax = (10-%) amp./cm?. (5) 
d* 


Equation (5) predicts that for a sample of unit 
length (d=1 cm) with 100 volts applied, a 
mobility of the charge carriers, b>=100 cm/sec. 
per volt/cm,* and an optical dielectric constant 
of the order of 10, a current density jmax~ 107° 
amp./cm? would be obtained. However, the 
relation between current density and voltage 
for Tl,S is given by Ohm’s law (Part I,! Section 
1B, Eq. (1)) rather than by Eq. (5), and a current 
density of about 1 amp./cm? would be observed 
under the conditions stated. 

A more likely assumption is that thermal 
ionization creates an equilibrium density of 
electrons and holes in the material : 


pa = poe @!*7, (6) 


where Q=ionization energy and pa=equilibrium 
density of electrons=equilibrium density of 
holes. If these charge carriers can traverse the 
compound with mobilities 6, and 6, without 


* This value corresponds to the maximum excess electron 
mobility in the alkali halides (A. Smakula, Goettinger 
Nachrichten Math. Phys. KI. 1, 62 (1934)), to the hole 
mobility in CusO (E. Engelhard, Ann. d. Physik 17, 501 
(1933)), and to the order of magnitude of the electron 
mobility in metals. 


space charge formation, a dark current density, 
ja, Would result: 


Ja= pat (be+dy). (7) 


The motion of charge carriers in semicon- 
ductors is normally an intermittent process.‘ 
The electron or hole travels a mean distance w in 
the field direction for an average time, t;, and 
is then trapped for a time, ¢;. The mobility of 
the carrier is thus given by 


b=Wo/(ts+ts), (8) 


where #)=W/E. The activation energy, u, re- 
quired to release the trapped carrier is supplied 
by thermal vibration; hence, ¢, will depend on 
temperature in accordance with Eq. (9): 


ty = toe"!*?, (9) 


For ¢,>>t;, which is usually the case, the mobility 
may be written 


b=boe“/*?, (10) 


It follows from (6), (7), and (10) that the current 
density can be expressed as 


ja=preE -e-2*? (bey exp [ —ue/kT] 
+bro exp [—u,/kT]). (11) 


(Subscripts e and h refer to electrons and holes, 
respectively.) 

It will be recalled that unoxidized T1.S (sample 
No. 2) proved to be an excess conductor and that 
on oxidation the resistance of the material 
increased rapidly, while the thermo-voltage 
changed its sign to indicate defect conductivity 
(Part I,' Section VI, Fig. 17). Consequently, 
pure thallous sulfide appears to be an ‘“‘intrinsic”’ 
semiconductor in which electrons and holes are 
mobile, the mobility of the former dominating 
(b.>by). 

In Part I,' Section V, Eq. (4), an empirical 
relation between conductance, Go, and tem- 
perature had been found: 


Go=ae-U/kT, 


Also, Table II had indicated that the activation 


‘A. von Hippel, J. App. Phys. 8, 815 (1937); Phys. Rev. 
8, 605 (1940). 








372 A. VON HIPPEL AND E. S. 


energy, U, varied appreciably from cell to cell 
and that no simple relation existed between U 
and the state of oxidation. Equation (11) clarifies 
this situation somewhat. One of the exponential 
terms involving the mobility normally dominates. 
In the single exponential function remaining, 
U is not a simple activation energy but the sum 
of an ionization energy Q and an activation 
energy u. It is likely that Q and u will be affected 
by oxidation. In addition, there is some evidence 
that in the very thin photo-cell layers a and U 
are functions of the layer thickness. Moreover, 
the insulating compound, Tl.SO:, formed in the 
oxidation, will influence U and a in accordance 
with the amount formed and its location. There- 
fore, the random variation in U from cell to 
cell is not surprising. 

For the unoxidized material according to (11) 
the conductivity can be written: 


o = poebeo exp [(O+u.)/kT |. (12) 


On a thin layer of pure Tl.S(~500A) values of 
¢=10-* ohm™ cm and of U=Q+4u,=0.54 elec- 
tron volt were obtained. The magnitude of po, 
which corresponds to the number of Tl atoms 
per cm? of thallous sulfide, is 1.110”. Hence 
the maximum mobility of the electrons at room 
temperature is be9=57 cm/sec. per volt/cm, in 
good agreement with values obtained from other 
systems.* 

From a consideration of the available space 
in the T1l,S structure,* it would appear that the 
most likely sites for the solute oxygen atoms 
are in spaces interleaving the TI-TI layers. This 
disposition provides twenty-seven possible loca- 
tions in the unit cell for accommodating oxygen 
atoms—the same as the number of sulfur atoms 
already present. Such strategic location of the 
oxygens permits them to exert a profound in- 
fluence on the conductivity. Owing to their large 
electron affinity, the oxygen atoms serve as 
electron traps, but the electrons apparently can 
be released again by thermal agitation. With 
increasing oxygen content the mobility of the 
excess electrons must rapidly decrease and the 
mobility of the holes, while somewhat impaired 
by lattice distortion, begins to dominate (Part 1, 
Section VI, Fig. 17). 


* The advice of F. G. Chesley has been most helpful in 
these structural considerations. 
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II. QUALITATIVE DISCUSSION OF THE 
PHOTO-EFFECT 

Two alternative explanations might be offered 
for the fact that oxygen photo-sensitizes thallous 
sulfide. A new, photoelectrically active, ,absorp- 
tion band may result from the oxidation, or the 
absorption of the Tl:S proper may lead to a 
photoelectric effect when oxygen is present. The 
transmission measurements (Part I,! Section \ 11) 
suggest the latter to be the case. 

Since the flash response characteristics show 
that a time of 510-5 sec. is required for the 
build-up of the photo-current (Part I,! Section 
VIIIB), it appears that initially an excited state 
is formed which persists for a short time before 
giving rise to longer-lived disturbance centers 
responsible for the secondary photo-effect. The 
details of this process are probably the following. 
By the absorption of a light quantum an electron 
is transferred instantaneously into an excited 
state. The creation of this state alters the poten- 
tial field of the neighboring atoms; the nuclei 
swing into new equilibrium positions, and the 
lifetime of the excited state is lengthened to 
about 10-5 sec. The electron can then either 
return to its original state accompanied by re- 
alignment of the lattice, or it can be ejected 
leaving behind a positive hole. 

In insulators the latter process would initiate 
only a primary photo-effect: that is, if electron 
and hole were to migrate through the material 
to the electrodes, a maximum _ photoelectric 
yield of one electron per ionizing light quantum 
would be recorded. In semiconductors. an addi- 
tional secondary photo-current should be ob- 
served, if the electron and hole do not leave the 
substance at the same rates. 

As is known from gas discharges, one positive 
ion moving slowly from anode to cathode can 
allow many fast electrons to pass in the opposite 
direction because it compensates at any moment 
the space charge density of the electron current. 
In the same way a slow hole moving from anode 
to cathode in pure T1.S (6.>),) would allow 
b./b, “extra” electrons to flow through the 
material. Unoxidized thallous sulfide has been 
found (Part I,! Section VIIIB) to exhibit such a 
photo-response, although the, yield obtained is 
lower than that of the oxidized material by 4 
factor of 104 to 10°. 
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The presence of an oxygen atom, with its 
large electron affinity, enormously increases the 
probability that a photoelectrically excited elec- 
tron will escape from its parent atom leaving 
behind a positive hole** and forming a relatively 
immobile O- ion. This ion represents the dis- 
turbance center producing the secondary photo- 
effect, and until it recombines with a positive 
hole, a number of mobile holes can traverse the 
material. The rate of disappearance of centers is 
proportional to the total number of O- ions and 
also to the total number of positive holes present 
at any time. 

The possibility of effecting conversions from 
one wave-length to another by means of the a 
factors in the plots of AJ vs. light intensity 
(Part I,! Section VD) and the lack of dependence 
of the shape of the decay characteristic on the 
wave-length of the monochromatic light flash 
(Part I,! Section VIIIB) indicate that the same 
kind of disturbance center is produced at differ- 
ent wave-lengths. 


Ill. QUANTITATIVE DISCUSSION OF THE 
PHOTO-EFFECT 


A. Response to Steady Illumination 


A calculation of the photo-response according 
to the picture just presented consists of an 
evaluation of the charge density of electrons and 
holes prevailing in the material and its change 
as a function of intensity of illumination. 

When thermal equilibrium is established in the 
dark, the rates of formation, dp,;/dt, and of 
recombination, dp,/dt, of the charge carriers are 
equal. Since the recombination rate is propor- 
tional to the number of electrons and to the 
number of holes present at any time, at equi- 
librium these rates will be proportional to the 
square of the charge density, py. 


dp; dp, 2 
—— = Bp2’. (13) 
dt dt 


(B=probability of recombination.) 

Under constant illumination a new balance 
between the rates of recombination and forma- 
tion is established and the charge density in- 


** The positive holes can probably be identified with 
TI* ions located at lattice points. 


creases from pg to p. Hence, 
Bp? = Bpa’+cJ, (14) 


where c represents the number of charge carrier 
pairs created per second per unit of volume 
under an illumination, J, of 1 wwatt/cm®. The 
increase in charge carrier density is therefore 


i 
2 


P 
Ap = p— pa= — pat (v7 + J) 
B 


and 


Ap cJ \3 
= -14(14 ). 
Pa Bp.” 


Since the photo-current density is given by 
Aj=ApeEbd,, (17) 
while the dark current density is given by 
ja=paeE(b.+b,), (18) 


the ratio of light current to dark current density 
follows from (16)-—(18): 


Aj bs cJ \3 
HL 2)] os 
Ja be +b; Bp.* 


In the sensitive, oxidized cell where b.<«b,, (19) 
reduces to 


Aj cJ \3 
= -14(14 ). (20) 
Ja Bp? 


. 


Values of Aj/ja calculated from (20) for a 
sensitive cell agree within experimental error 
with measured values of AI/Io, as can be seen 


in Table I. 


TABLE I. Agreement between measured and calculated 
response to continuous illumination. 


Illumination Al /Io Aj/ja 
Source » watts/cm? observed calculated 





White light 9.35 m 2.30 
23.0 3. 4.02 
42.1 y 5.72 
84.5 8.46 
214 14.1 
758 27.2 


Aao 


0.0877 
0.175 
0.263 
0.351 
0.438 
0.527 


2 


-165 
315 
440 
560 
.670 
.778 


ssssss FFye 
NOt 
ROO 








* The value of c/Bpa* at any other wave-length for this cell is given 
by (4.07) X(a factor). (Part I (reference 1), Section VD.) 
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PHOTORESPONSE , A j , MICROAMPERES/CM* 


—— EXPERIMENTAL 
———CALCULATED 


MINUS 


TIME, MILLISECONDS 


Fic. 1, Experimental and calculated decay characteristics 
(ja = microamperes/cm?). 


For incompletely oxidized cells (19) fits the 
data better than (20) because the excess con- 
ductivity is still appreciable. 


B. Response to Flash Illumination 


Under flash illumination the rise of the current 
occurs as soon as the charge carrier pairs are 
created. The number of these charge carriers 
should be proportional to the number of quanta 
absorbed. If the rate of formation of the pairs 
is much faster than the rate of their recombina- 
tion, the peak of the current transient should be 
proportional to the light intensity. 

The decay characteristic can be calculated 
from the differential equation expressing the 
rate of disappearance of charge 


dp/dt = B(pa?—p’). (21) 


Introducing Ap, the increase in charge carrier 
density above pa, 


dAp/dt = — BAp(Ap+2pa), (22) 
rewriting (22) 
1 1 
dap(—-—) =—2p.Bdt, (23) 
Ap Ap+2pa 


AND E. §S. 





RITTNER 
and integrating (23) 
Ap App 
2patAp 2patApy 





exp (—2paBt) 


(24) 


(Ap,»=maximum charge carrier pair density pro- 
duced by the light flash at t=0 which corre- 
sponds to about 5 10~-° sec. after occurrence of 
the flash), the equation for Ap is obtained: 


Apy*pa exp (— 2paBt) 
- 5 
2pat+Ap,(1—exp (—2paBt)) 


App exp (—2paB) 
= ~ (25) 


App 
1+—(1—exp (—2p,B?)) 
2pa 





A 





The photo-current density, Aj, plotted loga- 
rithmically as function of time, should therefore 
follow the law 


In Aj=In Aj,—2paBt 


App ; 
—In | 1+—(1—exp (—2p2Bt)) |, (26) 
2pa 


where Aj,=Ap,eEb,. 

Differentiating (26), the slope is found to be 
—(2patAp,)B at t=0 and —2p.B at t=~. 
Hence, the final slope of the decay curve is 
independent of light intensity while the magni- 
tude of the initial slope increases with the 
illumination level, in qualitative agreement with 
the experimental data (Part I,! Section VIIIB, 
Fig. 20). 

The constants needed to calculate the decay 
characteristic from (26) can be obtained in the 
following way. The quantity 2paB is equal to 
the final slope of the curve, and the extrapolation 
of this straight line to t=0 gives an intercept 


App App 
am aj / (1+), 

2pa 2pa 
from which Ap,/2pa can be computed. Substi- 
tuting the known value, Aj;, at t1(t2>1/2paB) in 

(26), the initial intercept, Aj», is obtained. 

Agreement between experimental decay char- 
acteristics and those calculated from (26) is not 
very good at the initial portion of the curves 
(Fig. 1). The difference between experimental 


(27) 





(24) 
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and theoretical characteristics is shown by the 
dotted curve in the figure, the height and shape 
of which suggest that an appreciable fraction, 
[Ap,’/(Appt+Ap,’) ], of the total number of dis- 
turbance centers decay at a faster rate than 
that predicted by (26). 

This situation may be accounted for by the 
following hypothesis. When the oxygen content 
of the solid solution is high, the electron affinity of 
an oxygen atom will be lowered by the presence 
of a close oxygen neighbor and the probability 
of recombination of the corresponding O- ion 
with a TI* hole will be increased to B’. In general, 
the number of oxygen close neighbors surround- 
ing an oxygen atom will vary statistically from 
zero to six and a similar statistical distribution 
of time constants is to be expected. It is the 
longest of these time constants, however, that 
exerts by far the profoundest effect on the 
response to steady illumination. 


C. Lifetime of the Photoelectric Centers 


On the basis of the preceding discussion, the 
average lifetime, r, of the disturbance centers 
can now be calculated. 

If an initial charge density, Ap», is created by 
a light flash, 7 is defined by the equation 


* Apdt 


T= ’ 
0 Ap» 


(28) 


or according to (25) 





“ 2pa exp (—2paBt)dt 
=f 2patAp,(1—exp (—2paBt)) 
By introducing x =exp (—2p,aBt) and 
dx = —2paB exp (—2paBt)dt, 


Eq. (29) can be rewritten 





T= 


- f dx 
B Jo 2pat+Ap,(i —x) 


1 
=— [In |2pat+App(1—x)} Jo? 
BApy 


1 Ap 
= In (1+). (30) 
BAp, 2pa 


The average lifetime, therefore, decreases with 


increasing light intensity. For very small light 
intensities (Ap,&K2pa)r becomes 


1 1 
T= =— (31) 
2paB K 


(for definition of x see Part I,! Section VIIIB). 


D. Frequency Response 


Runge and Sewig® have proposed a theory to 
account for the photo-effect in thallous sulfide 
photo-conductive cells based on the assumption 
that light creates disturbance centers which 
decay exponentially with time. This theory is 
inadequate to explain most of the data obtained 
in our work.' However, their assumption is 
valid at very low intensity levels, as is shown by 
flash response results (Part I,!} Section VIIIB, 
Fig. 20). Hence, the equation derived by these 
authors for the photo-current produced by a 
square-wave modulated light signal can be em- 
ployed to calculate the frequency response of the 
cell at the low intensity levels employed in the 
standard test set.*? In this equipment the d.c. 
component and all harmonics of the photo- 
current above the fundamental frequency are 
eliminated by the amplifier system. Therefore, 
the relationship between the instantaneous photo- 
current, Ai, and the time, ¢; the frequency, w/27; 
the d.c. current, Ai,, produced by steady illumi- 
nation of the peak value; and the time constant, 
1/x, is given by 


_ 2ndty w K 
At = ( sin wt+— cos ot) (32) 


7 w?+ x? w?+K? 





or 
2xAty 
At =—————- sin (wt+6), (33) 
1(w?+x?)3 


where tan @=x/w. The peak amplitude, A, of 
the current is thus given by 


A =2Khiy/a(w?+x?)}, (34) 
or the r.m.s. voltage, S, developed across the 


5]. Runge and R. Sewig, Zeits. f. Physik 62, 726 (1930). 
( s yi L. Hole, N.D.R.C. Report No. 16.4—Special No. 9 
1944). 
7™L. N. Holland, N.D.R.C. Report No. 16.4—Special 
No. 6 (1944). 
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Fic. 2. Relative response vs. frequency for different 
values of the time constant. 


load resistor, R,, of the input circuit amounts to 


2xAi,»R1 
S=— —--, 
m(w?+ «?)3(2)3 

The quantity «/(w?+x’)', representing the 
relative response with frequency, has_ been 
plotted in Fig. 2 as function of the modulating 
frequency for different values of the time con- 
stant, 1/x. For cells which decay quite slowly 
(x*<w?) the relative response is inversely pro- 
portional to the frequency and independent of x; 
while when «x becomes of the same order as w, 
the shape of the response characteristic is highly 
dependent upon the magnitude of the time 
constant. 

It follows from (35) that the frequency re- 
sponse of a cell can be calculated if x is deter- 
mined from flash response data and if the photo- 
current, Az,, resulting from the exposure to an 
unmodulated signal flux of one microhololumen 
is measured.* Tables II and III show that satis- 
factory agreement between calculated and meas- 
ured values of S is obtained. 


(35) 


* Ai, can also be calculated from Eq. (27) of Part I 
(reference 1) using the value of c/Bp,? obtained from 
measurements at higher levels of illumination and cor- 
recting for the fact that only a fraction of the cell surface 
is illuminated in the test set. 


AND E. S. 


RITTNER 


It is of interest to consider the change in cell 
response resulting from a change in the time 
constant, 1/x, keeping all other factors constant. 
Since Ai, is proportional to the average lifetime 
of the disturbance centers and hence propor- 
tional to 1/x at low levels of illumination, 
(w’+ x?) should be representative of the peak 
amplitude, A. Values of (w?+x?)- for various 
values of 1/« and for w/27=90 cycles are shown 
in Table IV. The output at 90 cycles increases 
with time constant; however, there is little to 
be gained in increasing 1/x beyond the values 
of the time constants of the majority of cells 
made by present methods. Improvement of the 
frequency response, which requires a reduction 
in the time constant, can be obtained only at 
the expense of the absolute signal level. 

Since x=2paB, any means which can _ be 
utilized to increase pa (or the apparent pa) can 
be employed in reducing the time constant. 
Such means include exposure of the cell to high 
temperatures or to background unmodulated 
illumination. 


E. Photoelectric Yield 
The space charge compensating action of the 


disturbance centers during their average life- 


TABLE II. Comparison of calculated and measured signal 
voltages at 90 cycles. 








Signal, db 
(reference level 
0.1in volt =0 db) 
Calcu- 
lated Measured 


Aip, RL, 
“amp. megohms 





50B 3. “3 0.20 . 101 98 
54B : 2 0.07 108 108 
84B 8. 0.02 117 115 
40B 5. a 0.025 117 114 
66B ‘ ‘ 0.015 111 110 
18A 3. * , _ 58 








** Aiy too small to measure; calculated value from measured signal 
=3.6 X10~4 wamp. 


TABLE III. Comparison of calculated and 
measured frequency response. 








Frequency, 
cycles /sec. 30 90 180 360 720 


S,db calculated 107 101 95 89 83 
measured 104 98 93 86 «81 





54B S,dbcalculated 114 108 102 96 90 


measured 112 108 103 94 88 


S, db calculated* 58 58 57.5 56 
measured 58 58 58 56 


18A 








* Calculations made relative to measured response at 90 cycles. 
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time, 7, allows z holes to traverse the length, d, 
of the semiconductor where 


2=rb,E/d. (36) 


If the fraction of the incident radiation absorbed 
by the layer is f, and if the photoelectrically 
excited electrons escape from their parent atoms 
with a probability, P, the quantum yield is 
given by 
n=fPrb,E/d. 
The photoelectric yield is a function of tem- 
perature for several reasons. As is shown by the 
data of Table IV, Part I,' the variation of the 
peak of the flash response curve with temperature 
can be expressed by an equation of the form 


AG,/Go=p—aT, (38) 


where p=2.0 and g=.0222 for T in °C and for 
the particular cell measured ; or 


AG, = pGo—qGoT 


(37) 


(39) 


(See (4), Part I.!) Hence, AG, increases at first 
and later decreases with temperature, becoming 
zero at 90°C. Since AG, is proportional to fP, 
this product also varies with temperature in 
accordance with (39). ' 

The density of charge carriers, py, increases 
rapidly with temperature according to (6) ; there- 
fore, if B is constant, the time constant for low 
light intensities should increase exponentially 
with the reciprocal of the absolute temperature. 
This has been experimentally verified (Part 1,' 
Fig. 22); however, the slope of the plot of 
log 1/x vs. 1/T was steeper than the corre- 
sponding plot of log Go vs. 1/T. This suggests 
that B is not a constant but varies with tem- 
perature in the following way: 


B= Bye!*7, (40) 


Finally, }, should increase somewhat with 
temperature (10), but the activation energy, un, 
is so small (see IJIF) that this effect is of minor 
importance. 

The 6 factors (Part I,' Fig. 15) representing 
the over-all yield as a function of temperature 
are defined by the relation 


(for \=constant) 


b=nrex/n300°K 


=fpxKP pext reK/fs00°KP300°K7300°K: (41) 
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TABLE IV. Dependence of (w?+x«?)~4 on 1/«x at 90 cycles. 








(w? +-42)-4 


4x 107° 

3.9X 10-4 
8.7 X 10-3 1X 10-3 
1.6X 10-3 4x10 
1.76 X 10-8 4x10 
1 


1/x, sec. 


4x10-5 
4x 10-4 





.167 X 10-3 90 











F. Calculation of Fundamental Constants 


If the value of one constant is assumed, it is 
possible to calculate the values of the other 
fundamental constants from experimental meas- 
urements. This has been done for the sensitive 
cell No. 84B. 

The hole mobility, },, can be obtained from (37) 
using measured values of quantum yield (y=1.8 
electrons/quantum at \=0.6u and E=124 volts/ 
cm), time constant at low levels of illumination 
(7 =9.13 X 10~* second, from final slope of Fig. 1), 
distance between alternate grids (d=0.19 cm), 
fraction of the incident radiation absorbed by the 
layer (f=0.75 at 0.64 from absorption measure- 
ments), and assuming the probability that a 
photoelectrically excited electron will escape from 
its parent atom, P, to be about 0.4: 

cm _ ,volt 
b, =0.96— / —-. 
sec. cm 

The maximum charge carrier pair density, 
Ap», produced by a single light flash (of unknown 
absolute intensity) can be computed from the 
calculated initial intercept of the flash response 
curve (Aj,=347 X10~* amp./cm? at \=0.9y and 
E=215 volts/cm), the calculated value of dy, 
and (26). 

Ap,»=1.0X 10" cm~*. 


The extrapolation of the straight line portion 
of the flash curve yields a value of Ap,/2pa=0.22. 
Hence, 

pa=2.4X 10" cm~*. 


Since the final slope, 2paB, is 109.5 sec.—', the 
recombination probability, B, is 
B=2.3X10-" cm? sec.~. 


Similarly, the average probability of recom- 
bination, B’, of the faster decaying centers is 
about 

B’~1.6X10-" cm? sec.—, 
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and these centers comprise about 43 percent of 
the total number. (See initial intercepts of 
Fig. 1.) 

From the calculated value of pg, the known 
values of po (1.110”) and T (300°K), and (6), 
the thermal ionization energy can be calculated. 


Q=0.51 electron volt. 


This value agrees with the measured thermal 
activation energy, U, within experimental error. 
Hence, 


u,= U—Q=0. 


The quantum yield measured under steady 
illumination yields a value of c/Bp,=11.0, 
which together with the values of B and pa gives 
the number of charge carrier pairs created per 
second per unit of volume under an illumination 
of 1 wpwatt/cm? (A=0.6yz). 


c=1.4X10'* cm sec. pwatt cm?. 


This corresponds to an equilibrium charge den- 
sity created by this illumination of Ap=5.9 
X 10! cm~. 

The constant, c, can also be calculated di- 
rectly from the number of effective quanta 
absorbed in the photo-sensitive layer of thick- 
ness 6 (6=5X10-5 cm). 


5.03 X10"A,fP 
C= 
6 





= 1.8 10'* cm™ sec.—! pwatt-! cm*. (42) 
In view of the uncertainty in the magnitude of P, 
this result is in good agreement with the previous 
calculation. 

Finally, the capture cross section, A, of an 
oxygen ion defined by 


A=B/v (43) 


can be obtained if an estimate is made of the 
“instantaneous” velocity, v, of a thallium hole. 
The maximum frequency with which a TI* hole 
can be exchanged from thallium to thallium 
neighbor by thermal vibration can be calculated 
to a first approximation from the Debye tem- 
perature of thallium metal: 


Vmax = kO/h 


= 2.110" sec.—. - 
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where k=Boltzmann constant, 4=Planck con- 
stant, and @= Debye temperature = 100°K. Mul- 
tiplying ymax by the average thallium —thallium 
distance in TlS (410-8 cm), the maximum 
velocity, Umax, of a thallium hole is found to be 


Umax = 8.4 X 104 cm sec. 
=v for u4,=0 


yielding a capture cross section of 
A=2.74X10-" cm’, 


an area which is unreasonably small. This indi- 
cates that many thallium holes collide with 
oxygen ions without recombining. Evidence has 
already been cited to support this supposition 
(see (40)), an activation energy, s=0.11 electron 
volt, having been measured on a different cell 
(Part I,! Section VIIIB). If this value of s also 
applies to cell No. 84B, the maximum probability 
of recombination is 


By=1.61X10-" cm? sec.“, 
leading to a value of 
A=19.110-* cm?, 


which corresponds to a circular cross section of 
radius 2.46A. This is a reasonable order of 
magnitude as the distance from center to center 
of a TI* ion and O- ion in contact is about 2.5A. 
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Notes on the Infra-Red Spectrum and Molecular Structure of Ozone 
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The nature of the fine structure of the band at 14.2y in the absorption spectrum of ozone 
has been determined with prism and grating spectrometers. These data, together with the 
fine structure of the 9.574 band (previously resolved in the solar spectrum) render possible an 
unambiguous choice between the divergent views currently held with regard to the structure 
of the ozone molecule. It is concluded that the form of the isoceles triangle is acute, with the 


apex angle in the neighborhood of 34°. 


I, INTRODUCTION 


UMEROUS attempts have been made to 
analyze the infra-red spectrum and de- 
termine the molecular structure of ozone. 

Prior to the definitive (although low disper- 
sion) experiments of Hettner, Pohlman, and 
Schumacher! in 1934, these attempts were un- 
successful, mainly because 


(1) they were made in ignorance of the existence 
of the ozone band 4t 14.2, and 

(2) they mistakenly identified as due to ozone 
the intense bands of nitrogen-pentoxide 
(N2O;), a contamination difficult to avoid in 
the preparation of ozone. 


Since 1934 the attempts at solving the ozone 
problem have been inconclusive because of the 
persisting lack of accurate information regarding 
the band envelopes and fine structures. 

The history of the problem prior to the work 
of Hettner, Pohlman, and Schumacher’ is given 
by Gerhard,? and since that time by Simpson.* 
It is sufficient for our present purpose to point 
out that Hettner, Pohlman, and Schumacher! as- 
sumed a Q-branch for the 9.574 band. The en- 
velopes of the bands at 4.754 and 14.24 were 
plainly Q-branch and doublet, respectively. On 
the assumptions of central forces and the fre- 
quency assignment 71(4.75y), ve(9.57u), v3(14.2) 
this led to an apical angle of 39 degrees. 

Low dispersion solar spectra show, however, 
an apparently doublet envelope for the band at 





_'G. Hettner, R. Pohlman, and H. J. Schumacher, 
Zeits. f. Physik 91, 372 (1934). 

*S. L. Gerhard, Phys. Rev. 42, 622 (1932). 

*D. M. Simpson, Trans. Faraday Soc. 41, 209 (1945). 


9.574. Adel, Slipher, and Fouts‘ showed that, 
using this identification, imaginary values for the 
apical angle are obtained on either central or 
valence forces, provided the bands at 4.75y, 
9.574, and 14.2 are still* regarded as funda- 
mentals. 

Continuing the assumption of doublet struc- 
ture for the 9.574 band (which we believe to be 
a mistaken one), but substituting the weak band 
at 5.75u for the strong one at 4.75u, and giving 
preference to a valence force field, Penney and 
Sutherland® were led by the frequency assign- 
ment 71(9.57u), v2(14.2u), v3(5.754) to an apical 
angle of 164 degrees. 

In a critical discussion of the electronic struc- 
ture of the ozone molecule Mulliken® favors the 
acute angle form of the molecule designated by 
Hettner, Pohlman, and Schumacher.' But re- 
cently Shand and Spurr’ have concluded from 
electron diffraction experiments that the apical 
angle of the ozone molecule must be obtuse, 
127+3 degrees. 

The correctness of this last result was assumed 
by Miss D. M. Simpson.* Continuing the tradi- 
tion of the doublet envelope for 9.574 she finds 
the most probable assignment of frequencies to 
be v1(5.75u), vo(9.57u), v3(4.75y). 

In deciding between an acute and an obtuse 
form for the ozone molecule, the envelope form 
of the 9.574 band is crucial. And since, in all 
analytical attempts to date, this form is at least 
somewhat uncertain, it is desirable to obtain 





4A. Adel, V. M. Slipher, and O. Fouts, Phys. Rev. 49, 
288 (1936). 
5 W. G. Penney and G. B. B. M. Sutherland, Proc. Roy. 
Soc. A156, 654, 678 (1936). 
6 R. S. Mulliken, Rev. Mod. Phys. 14, 204 (1942). 
( aame Shand and R. A. Spurr, J. Am. Chem. Soc. 65, 179 
1943). 
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new and independent evidence. This has been 
achieved by comparing the fine structure of the 
doublet band at 14.2 with the fine structure of 
the band at 9.57y. 

The fine structure of 9.57y is already at hand 
in the infra-red solar spectrum.® It remained to 
resolve the band at 14.2uy. 


II. EXPERIMENTAL 


Ozone was prepared according to the method 
of Goldstein,® and the 14.24 band was examined 
both with a NaCl prism and a 2400-line grating. 
The prism curve, Fig. 1, is remarkably suggestive 
of the perpendicular band of a symmetrical top 
molecule.* Asymmetry in the lines of Fig. 1 sug- 
gests a complex line structure, and this is con- 
firmed by the grating analysis, Fig. 2. 

The structures of the bands at 9.574 and 14.2u 
are fundamentally dissimilar. And since 14.2y is 
a near-perpendicular band, 9.574 must be a 
near-parallel band (Q-branch). 

As a further aid in discussing the ozone spec- 
trum the envelopes of the bands at 3.28u and 
3.57u have been re-examined, Fig. 3. 


III. DISCUSSION 


The first and the most crucial point in the 
interpretation of the infra-red spectrum of ozone 


8 A. Adel, Astrophys. J. 94, 451 (1941). 

® E. Goldstein, Ber. d. d. Chem. Ges. 36, 3042 (1903). 

*We are indebted to Mr. E. A. Boettner of the 
Wyandotte Chemicals Corporation for the use of their 
prism spectrometer in obtaining the data for this curve. 
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is the identification of the fundamental bands. 
The investigation of Hettner, Pohlman, and 
Schumacher! showed that the most intense re- 
gions of absorption are those at 4.75, 9.57, and 
14.24. From 16 to 27u (the limit of their observa- 
tions) ozone is almost completely transparent 
and it is unlikely that any fundamental bands 
would lie beyond 27. It is therefore natural to 
identify the three regions just mentioned with 
the three fundamental frequencies of the mole- 
cule. An alternate scheme which has been con- 
sidered by some investigators would be to take 
the band at 9.57y as a fundamental but to iden- 
tify the 4.75y region as its first overtone. The 
second overtone must then be assigned to the 
weak band at 3.28u. The difficulty with this 
scheme is that the three frequencies, namely, 
1043, 2108, and 3050 cm— cannot be placed in a 
series of the type demanded by the theory of 
anharmonic oscillators. This is immediately evi- 
dent since 2108 is greater than twice 1043, while 
3050 is less than three times 1043. 

For these reasons it seems very probable that 
the three fundamental frequencies of ozone are 
v;=2108, ve=1043, and v3=705 cm™, corre- 
sponding to the three most intense regions of 
absorption. 

The second crucial point is concerned with the 
fine structure of the bands. The new experi- 
mental evidence reported in this paper sheds a 
revealing light upon the problem. The long wave- 
length band at 14.2 is found to have a structure 
consisting of widely spaced lines, each of which 
is multiple. The average value of the coarse 
spacing is about 5.7 cm~ and is fairly regular 
although there is a rather large convergence. 
The magnitude of the coarse spacing indicates 
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unequivocally that the electric moment corre- 
sponding to this frequency must oscillate along 
the middle axis of inertia. (The oscillating elec- 
tric moment of a triatomic molecule can only 
be directed along either the least or the middle 
axis of inertia. When it is directed along the 
middle axis, the rotational selection rules predict 
a widely spaced series of lines characteristic of 
the rapid rotations about the least axis of inertia. 
When the electric moment oscillates along the 
least axis of inertia, these transitions are 
forbidden.) 

The band at 9.574 has been examined by 
Adel® in the solar spectrum. Its structure was 
found to consist of fairly regularly spaced lines 
with an average spacing constant of around 1 
cm~!. Gerhard’s? measurements on the 4.75y 
band likewise reveal a very narrow spacing which 
is also of the order of 1 cm~!. Thus for both of 
these bands the electric moment must oscillate 
along the least axis of inertia. 

These facts admit of only one interpretation. 
The apical angle of the ozone molecule must be 
less than 60 degrees. This conclusion is the re- 
sult of broad and general considerations. A 
further examination serves to confirm. it in detail 
and leads to the following points. 

1. The spacing of the fine structure lines in a 
band where the electric moment oscillates along 
the least axis of inertia is of the order of Av 
=(h/8m’c)(1/Ia+1/Ip). (Of course this mole- 
cule is an asymmetric rotator and the spacing of 
the fine structure lines must be irregular. How- 
ever, when two of the moments of inertia, in 
this case the greatest and middle moments, are 
approximately equal, the above formula gives a 
rough representation of the fine structure.) Set- 
ting in the observed spacing of 1 cm=, one ob- 
tains 224Ip/(I4 +I) 56 X 10-* g cm’. 

The spacing of the fine structure lines in a 
band with changing electric moment along the 
middle axis of inertia is of the order of Av=(h/ 
82°c)(2/Ie—1/I4—1/Ig). This approximation is 
subject to the same limitations as those govern- 
ing the spacing in the other type of band. In 
the ozone molecule, the identity of the oxygen 
nuclei and the fact that their nuclear spin is 
zero means that every other line will be missing 
and the observed spacing will have twice the 
value given above. Setting in the experimental 
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value for Av, combining with the earlier result, and 
remembering that for a plane molecule J,4= 
Ip+TIc, one obtains [4264 10-*°, 73»—=50 10-”, 
and Ic=14X10-* g cm*. The corresponding api- 
cal angle is 34 degrees. It must be emphasized that 
these numerical results are extremely rough both 
from the experimental and the theoretical sides. 
Our estimate of their order of accuracy is +20 
percent. In order to obtain the moments of in- 
ertia more exactly the most promising procedure 
appears to be to resolve the complete fine struc- 
ture of the 14.24 band with a grating spectrom- 
eter. The lines can then be correlated with the 
energy levels of an asymmetric rotator and the 
moments of inertia determined. We hope to 
make this investigation in the near future. 

2. If the molecule of ozone has the form of an 
isosceles triangle with an apical angle of less than 
60 degrees, it is quite probable that the forces 
can be satisfactorily approximated by central 
forces rather than by valence forces. A calcula- 
tion of the force constants yields results which 
are substantially in agreement with those ob- 
tained by Hettner, Pohlman, and Schumacher. 
The apical angle turns out to be 39 degrees, the 
force constant between the oxygens forming the 
base of the triangle, K = 34X10° dynes/cm and 
the constant between the apical and one of the 
base oxygens, K’ =6.2 10°. 

3. The weaker bands of ozone may be ex- 
pressed as overtones and harmonics of the three 
fundamentals. These have been assembled in 

















382 A. ADEL AND D. 
TABLE I. 
Calc. Obs. 
Identification frequency frequency Obs. type 

V3 ——D 705 Doublet 

ve —Z 1043 Zero branch¢ 
vot vs 1748D 1740 Doublet® 

vy —Z 2108 Zero branch¢ 
vitvs 2813D 2800 Zero branch?4 
vitve 3151Z 3050 Zero branch?¢ 





« This band has often been thought to be of the doublet type. A study 
of Adel’s curve shows an intense region of absorption at 1041 cm™ 
which is undoubtedly a zero branch displaced through convergence. 

+’ While the resolution is not great, the curve of Hettner, Pohlman, 
and Schumacher distinctly favors a doublet structure for this band. 

¢ Gerhard’s work, as well as ours, shows that this is clearly a zero- 


branch band. : 
4 Figure 3 shows envelopes which more nearly resemble zero branches 
than doublets. The evidence is inconclusive for both of these overtone 


bands however. 


Table I. The first column gives the designation 
of the band and the second the calculated fre- 
quency and the predicted type (D=doublet 
band, Z=zero-branch band). The third column 
lists the observed frequency while the fourth 
describes its type. 


IV. CONCLUSION 


The present interpretation of the ozone spec- 
trum is both simple and direct and is, in almost 
all respects, self-consistent. The three most in- 
tense regions of absorption are identified with 
the three fundamental frequencies. The weaker 
regions are assigned to those overtone combina- 
tions which are normally the most intense, 
namely, ve+v3, vitvs and vi+v2. The calculated 
frequencies of the overtones are, in each case, 
slightly larger than the observed frequencies. 
This again is the usual situation. 

The fine structure of the bands leads to an 
apical angle of 34 degrees which is in good agree- 
ment with the apical angle of 39 degrees cal- 
culated independently from the positions of the 
fundamental bands. 

The dimensions of the molecule as roughly 
estimated from the fine structure are reasonable. 
The distance between the oxygens forming the 
base of the isosceles triangle is 1.0X10-* cm 
whereas the distance between the atoms in O2 
is 1.20810-°. The envelopes of the bands—a 
distinct doublet for v3 and zero-branch bands for 
vy, and ve—are in agreement with the predictions 
from the theory of the asymmetric rotator. 

The force constants are not unreasonable. 
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The force between the base oxygens is about 5.5 
times that between a base and the apical oxygen. 
The force constant between the base oxygens 
would predict a vibration frequency of these 
two atoms taken by themselves of about 2700 
cm! whereas the vibration frequency of Oz is 
1580 cm™. It is well known that normal vibra- 
tions and interatomic distances follow the ap- 
proximate law r."w.=constant where n lies be- 
tween 2 and 3. In our case the figures lead to a 
value for 1 of 2.8. (The above argument must be 
regarded merely as an indication since the struc- 
tures of O2 and of O; may be quite dissimilar.) 

The Raman spectrum of ozone is very weak 
and shows no lines which can be identified with 
certainty.1*" A molecule with an acute apical 
angle would not be expected to possess a strong 
Raman spectrum whereas one with an obtuse 
apical angle would. 

There are only two experimental facts which 
do not appear to be in concordance with our in- 
terpretation of the spectrum. (a) The dipole 
moment of ozone is known to be small,!” 0.49 
x<10-'8. In general a molecule with an acute 
apical angle should exhibit marked polar quali- 
ties although in any particular case other cir- 
cumstances might tend to reduce the dipole 
moment. (b) The electron diffraction measure- 
ments of Shand and Spurr indicate an apical 
angle of 127 degrees. 

The experiment appears to have been care- 
fully performed and we are unable to suggest 
any explanation for the discrepancy between the 
electron diffraction and the infra-red results. It 
is certainly true, however, that the infra-red 
spectrum of ozone cannot be identified as be- 
longing to a molecule with an obtuse apical 
angle without the introduction of what seems to 
be an artificial series of assumptions and inter- 
pretations. It is believed that further precision 
measurements of the fine structure of the band 
at 14.2 will definitely establish the moments of 
inertia of the molecule. 


10 G. B. B. M. Sutherland, Proc. Roy. Soc. A141, 515 
(1933). 

1G. B. B. M. Sutherland and S. L. Gerhard, Nature 
130, 241 (1932). 

2G. L. Lewis and C. P. Smyth, J. Am. Chem. Soc. 61, 
3063 (1939). 
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Studies on Glass. 


XVIII. The Heats of Solution of Crystalline and Glassy Glucose; 
the Heat of Mutarotation of a-Glucose 


RoBERT D. ROWE AND GEORGE S. PARKS 
Department of Chemistry, Stanford University, California 
(Received March 12, 1946) 


A calorimeter has been developed for the rapid solution of a powdered solid in a suitable 
solvent and the accurate determination of the heat of solution during the process. This calo- 
rimeter has now been used to determine the heats of solution of two preparations of crystalline 
a-d-glucose in water, both with and without mutarotation, and also the heats of solution of 
two preparations of glucose glass. Relatively accurate thermodynamic data for the crystal-to- 
glass transformation have been computed from these solution heats. Also small, but significant, 
differences per mole in enthalpy (AH293;=102+18 cal.), in entropy (AS293=0.24+0.06 cal./ 
deg.), and in free energy (AF29s°=31+10 cal.) have been found for these two glass prepara- 
tions. A value of AH29g=—175 (+5) calories has been obtained for the mutarotation of a 
mole of a-glucose, dissolved in water, to yield the equilibrium mixture of the alpha- and 


beta-forms. 





INTRODUCTION 


HE measurements here described were car- 

ried out during 1937-38 in order to provide 

as accurate data as possible concerning the rela- 

tive heat contents of crystalline a-glucose and 

glassy glucose. The present results supersede the 

earlier, essentially preliminary, values of Parks, 

Snyder, and Cattoir.! Consequently, their ther- 

modynamic calculations for the crystal-to-glass 

transformation are now revised in the light of 
these newer data. 

When this investigation was started, it was 
also considered likely that the heat of solution 
of a sample of glucose glass would show varia- 
tions with the history of its annealing, i.e., the 
time and temperature of storage, prior to the 
solution process itself, as previous measurements 
of the optical rotatory power? and the viscosity* 
of this material had disclosed rather appreciable 
changes with aging. No such effect was dis- 
cernible, however, in the present study, probably 
because of the character of these solution experi- 
ments and the time requirements here involved. 

Incidentally, the data for the heat of solution 
of the crystalline a-glucose in water, both with 

'G. S. Parks, L. J. Snyder, and F. R. Cattoir, J. Chem. 
Phys. 2, 595 (1934). 

*C. J. Jacobs and G. S. Parks, an unpublished study at 


Stanford University. 


(1939) S. Parks and J. D. Reagh, J. Chem. Phys. 5, 364 


and without mutarotation, provide a new and 
relatively accurate value for the heat of muta- 
rotation. 


METHOD AND APPARATUS 


After considerable preliminary experimenta- 
tion, the “‘ordinary”’ method of calorimetry was 
adopted for measuring the heats of solution of 
the various glucose samples in water at approxi- 
mately 25°C. With our final apparatus the 
actual solution process of the glucose, properly 
powdered, took place in about twenty seconds 
and complete thermal equilibrium was then 
attained two or three minutes later. The cor- 
rected temperature changes were readily com- 
puted by the Regnault-Pfaundler method from 
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384 R. D. ROWE 
two series of temperature readings, taken for 
periods of twenty to thirty minutes prior to and 
immediately following each solution process. 

Our calorimetric assembly, as finally developed, 
is represented schematically in Fig. 1. The 
calorimeter proper, designated by the lettering 
AAAA, was a cylindrical copper container of 
about 450-ml capacity, which had been given a 
number of coats of Bakelite varnish both inside 
and outside to protect the metal from alkaline 
reagents. It had a specially shaped, streamlined 
bottom and a tightly fitting lid. The stirrer was 
a four-bladed propeller of monel metal, operating 
in a small copper cylinder (BBBB) immersed in 
the solvent liquid and set off center in such a 
way as to circulate all of the liquid around the 
calorimeter thermometer (7;) and the walls of 
the container. Evaporation was -prevented by 
the tightly fitting lid, the special grease-filled 
seal on the stirrer shaft, and a rubber gasket 
about the thermometer. This container, the 
calorimeter, was surrounded by a silver shield 
(CCCC), in the form of a cylindrical can with a 
removable lid, to reduce radiation and convection 
effects and thus insure a small, uniform rate of 
heat transfer to or from the jacket, which in 
turn was a large copper can (DDDD) with a 
removable top. Six Bakelite pegs (only three are 
indicated in Fig. 1) supported both the con- 
vection shield and the calorimeter firmly, but in 
such a way that the latter could be easily re- 
moved. A’ special doorway, not shown in this 
figure, was provided through all three lids so 
that samples could be added with a minimum of 
heat transfer. These holes were closed by a 
removable, hollow outside door which supported 
covers for the convection shield and the inner 
calorimeter lid. A small spring on the support 
held the cover of this hole in the calorimeter lid 
tightly in place so as to prevent evaporation of 
the solvent. 

The copper jacket, with the silver shield and 
calorimeter therein, was finally suspended within 
a fifteen-liter crock, which was filled with water 
at a temperature of about 24°C. A small rotary 
pump circulated the water in this crock, and both 
the jacket lid and hollow door were also filled 
with water. The entire apparatus was then 
placed in a thick-walled Celotex box with a 
removable Celotex cover. A 1/15 hp synchronous 
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motor, connected by means of a system of belts 
and pulleys, drove the calorimeter stirrer and 
jacket circulating pump at constant speed, the 
former at 505 r.p.m. and the latter at 600 r.p.m. 

The temperature changes within the calorim- 
eter were measured to 0.0005° with a very fine 
Beckmann thermometer, manufactured specially 
for us by Max Tischer and Company (I]menau, 
Germany) and calibrated at 0.2° intervals by 
the Physikalisch-Technische Reichsanstalt. The 
jacket temperature was observed with an ordi- 
nary Beckmann thermometer (72). These ther- 
mometers were set at proper temperatures against 
a 0°-50° thermometer which had been calibrated 
to 0.01° by the National Bureau of Stand- 
ards. This latter thermometer was also used to 
obtain the temperature of the samples and of 
the room during the experimental runs. Lens 
magnifiers were used in reading all thermometers. 

The calibration of the calorimeter in the 
22°-25° range was effectuated by a method of 
mixtures. In this calibration procedure a weighed 
sample of distilled water, initially at about 20°, 
was added to the calorimeter, which had been 
previously prepared by a partial filling with 
water and a temperature adjustment to 24° or 
25°. By means of a copper-constantan thermo- 
couple in conjunction with a White potenti- 
ometer, the temperature of this added water 
was measured to within 0.001° as it entered the 
calorimeter proper from a pipette. The final 
temperature of the calorimeter was then deter- 
mined on its own Beckmann thermometer, 
against which the thermocouple had been stand- 
ardized. Six such calibration experiments yielded 
a mean result of 62.0 calories per degree with a 
mean deviation of only +0.3 calorie. This result 
checked roughly with a prior estimate of 59.2 
calories per degree, exclusive of the heat capacity 
of the Beckmann thermometer, which had been 
made from the weights and known specific heats 
of the materials used in the construction of the 
calorimeter. 

MATERIALS 
Crystalline Glucose 

Two preparations of crystalline a-d-glucose 
were employed in the present study. 

The first material was a C. P. Pfanstiehl 
product, which was dried in a vacuum desiccator 
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over P.O; for six weeks before use. It was then 
estimated to contain around 0.01 percent mois- 
ture in view of Jackson’s observations.* An 
analysis by the polarimeter, based on the rotation 
data of Isbell and Pigman,> indicated that it 
was about 96 percent alpha- and 4 percent beta- 
glucose. This dried material was used in our 
initial measurements of the heat of solution of 
crystalline glucose and also in preparing the first 
sample of our glucose glass. Some of the undried 
crystals were likewise used in preparing the 
second sample of glass. . 

Our second crystalline material was some of 
Standard Sample No. 41, “Dextrose,” as pre- 
pared by the National Bureau of Standards. 
It was a fine powder and was guaranteed to 
contain not over 0.01 percent moisture. Its 
initial rotation indicated that it was 100 percent 
a-glucose within the limits of error of our 
polarimetric measurements. This material was 
used in our final measurements of the heat of 
solution of crystalline a-d-glucose and its heat of 
mutarotation. 


Glassy Glucose 


As in other studies,* the preparations of glassy 
glucose were obtained by heating the crystalline 
Pfanstiehl glucose in a flask in a paraffin bath. 
The crystals melted at approximately 146°C and 
were then maintained under a vacuum at 155- 
160° for about five minutes, in order to eliminate 
bubbles of air or water vapor. When the liquid 
had become fairly clear, air was admitted and 
the material was allowed to cool to the glassy 
condition. To facilitate use in the solution experi- 
ments, the glass so prepared was next reduced 
to a powder by passing it through an old- 
fashioned, home meat grinder. This powder was 
then stored in vacuum desiccators over P.O; 
until used. 

Two batches of glassy glucose were prepared 
for the present measurements. Glass I was made 
from the remaining Pfanstiehl crystals which 
had been previously dried over P.O; for six 


‘R. F. Jackson, Bureau of Standards Scientific Paper 
No. 293 (1916). 

°H. S. Isbell and W. W. Pigman, J. Research Nat. 
Bur. Stand. 18, 141 (1937). 

6G. S. Parks and W. A. Gilkey, J. Phys. Chem. 33, 1430 
(1929); G. S. Parks, L. E. Barton, M. E. Spaght, and 
J. W. Richardson, Physics 5, 193 (1934). 
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weeks. Glass I] was made from similar undried 
material, which, however, should become com- 
pletely freed of its water when heated above 
150°C. Half of this latter glass preparation was 
kept at room temperature; the remainder was 
kept in the laboratory ‘‘cold room”’ at about 3°C. 


THE HEATS OF SOLUTION 


A uniform procedure was followed in measuring | 
the heats of solution of these several glucose 
samples in water. Exactly 250 g of distilled water 
were first placed in the calorimeter and sub- 
sequently a portion of powdered glucose of 
45.0 (+0.2) g mass was dissolved therein. Thus 
the resulting solution was always 1.00 molal. 

The initial temperature of the calorimeter and 
solvent water in these experiments was about 
24°C. The solution process for the crystals, 
being distinctly endothermic, resulted then in a 
temperature drop of 1.8 to 2.0°. On the other 
hand, the solution process of the glass was some- 
what exothermic and yielded temperature rises 
of 0.8 to 1.0°. In the computations of the corre- 
sponding heat effects, the necessary values for 
the specific heat of water were taken from the 
International Critical Tables’; those for the 
glucose crystals and glass were derived from 
the preceding studies of Parks and co-workers.*® 
In order to put all results on a truly comparable 
basis, the value for the changes in enthalpy (AH) 
in the solution process were finally adjusted to 
the standard temperature of 25° (i.e., 298°K) by 
taking AC, =0.20 cal./g for the solution of the 
crystals and AC, =0 for the solution of the glass. 
These AC, values were based on the plausible 
assumption that the glucose in water constituted 
an almost ideal solution. The former value has 
also received excellent corroboration, subsequent 
to the completion of our work, from the measure- 
ments of Sturtevant,’ who has shown that for 
the solution of alpha-crystals to a final concen- 
tration of approximately 0.2 molal AC, =0.200 
cal./g between 25° and 35°C. Our second value 
of AC, =0 for the solution of the glass is of less 
importance, as the temperature adjustment to 


7 International Critical Tables (McGraw-Hill Book Com- 
pany, Inc., New York), Vol. VII, p. 232. 

8G. S. Parks, H. M. Huffman, and F. R. Cattoir, J. Phys. 
Chem. 32, 1366 (1928); G. S. Parks and S. B. Thomas, 
J. Am. Chem. Soc. 56, 1423 (1934). 

9J. M. Sturtevant, J. Phys. Chem. 45, 127 (1941). 
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TABLE I. Heat of solution of Pfanstiehl d-glucose 
crystals at 25°C. 








AH 2s for solution 





Experiment No. Solvent (cal. /g) 
P-1 Pure water 14.63 
P-2 Pure water 14.63 
Mean without mutarotation Pure water 14.63 (+0.02) 
P-3 0.025N HCl 14.64 
P-4 0.025N HCl 14.63 
Mean 0.025N HCl 14.63 (+0.02) 
P-5 0.0025N NHsOH = 13.99 
P-6 0.0025N NHsOH $13.97 
P-7 0.0025N NHs,OH = 13.96 
P-8 0.0025N NHsOH = 13.91 
Mean, incomplete mutarotation 0.0025N NH,OH § 13.96 (+0.03) 
P-9 0.0050N NH,OH 13.70 
P-10 0.010N NH.sOH 13.74 
P-11 0.020N NH,sOH 13.67 
Mean with mutarotation 13.70 (+0.03) 


Enthalpy change of mutarotation AH 293 = —0.93 (+0.04) 











25° in these cases involved an average interval of 
considerably less than one degree. 

The solution experiments with the Pfanstiehl 
crystals were essentially preliminary and explora- 
tory in character. These results are summarized 
in Table I. Experiments 1 and 2, carried out 
with pure water as solvent, yielded a value of 
AHg= 14.63 cal./g; and this result may be 
ascribed to the heat of solution of the material 
without mutarotation, as the crystals dissolved 
in less than twenty seconds and any subsequent, 
relatively slow mutarotation was automatically 
provided for in the derivation of the small, 
post-solution heating correction. The following 
pair of experiments with 0.025 normal HCl as 
solvent were expected to show at least some 
mutarotation effects, but apparently such a weak 
solution of hydrochloric acid is not an effective 
catalyst for the mutarotation process. Next, 
one to eight drops of concentrated ammonium 
hydroxide were added to the 250 g of solvent 
water to yield solutions ranging from 0.0025 up 
to 0.020 normal. The tabulated results show that 
in these cases the mutarotation process is ap- 
parently complete for our measurements with 
concentrations of, or above, 0.0050 normal. 

Following these exploratory experiments, the 
heat of solution of the high grade a-d-glucose, 
prepared at the National Bureau of Standards, 
was measured in a first pair of determinations 
with pure water as solvent and in a second pair 
with 0.0125N NH,OH as solvent. The results, 
representing the solution process both without 
and with mutarotation, are summarized in 


Table II. 


PARKS 


A comparison of these values for the heat of 
solution of a-glucose with most of the earlier 
data in the literature is fraught with uncer- 
tainty, as several of these earlier investigators 
either overlooked the mutarotation problem or, 
at least, made inadequate allowance for its 
effect. Thus Barry’s results,!® adjusted by us to 
25°C, are, respectively, AHo3= 15.04 (+0.09) for 
solution of one gram of glucose in pure water 
to yield 0.2 molal concentration and AH 2. 
= 14.90 (+0.16) for corresponding solution in 
1.64 normal hydrochloric acid, but he assumed 
that the heat of mutarotation was zero. Hendricks 
and co-workers" have also reported, without 
reference to possible mutarotation, 13.85 cal./g 
for the heat of solution of a-glucose in water at 
21.6° to yield a 0.017 molal concentration ; this 
result adjusts to AHo3= 14.53 cal./g. The result 
of Parks, Snyder, and Cattoir,! which corre- 
sponds to AHos = 14.72 cal./g, probably involved 
very incomplete mutarotation, although these 
workers recognized the possible importance of 
this effect. The best published work to the 
present date appears to be that of Sturtevant® 
who has reported AHo3= 14.22 (+0.02) cal. /g for 
the solution of a-glucose, presumably without the 
mutarotation effect, in water to yield an approxi- 
mately 0.2 molal concentration. His result is 
thus about 0.84 calorie below our present average 
for the non-mutarotated sugar. While his con- 
centration differs considerably from ours, it is 
really hard to explain such a large difference in 
heat of solution. 

Next, the heats of solution of preparations | 
and II of the glassy glucose were measured with 
0.0125N NH,OH as the solvent. This solvent 
was employed because the work of Jacobs and 


TABLE II. Heat of solution of pure a-d-glucose 
crystals at 25°C. 








AH 23s for solution 


Experiment No. Solvent (cal. /g) 
BS-1 Pure water 15.090 
BS-2 Pure water 15.041 
Mean without mutarotation Pure water 15.065 (+0.025) 
BS-3 0.0125N NHsOH 14.085 
BS-4 0.0125N NHsOH 1 


4.116 
Mean with mutarotation 0.0125N NHsOH 14.100 (+0.015) 


Enthalpy change of mutarotation AH 23 = —0.965 (+0.030) 








10 F. Barry, J. Am. Chem. Soc. 42, 1912 (1920). 
1 C, B. Hendricks, J. H. Dorsey, R. H. LeRoy, and 
A. G. Moseley, Jr., J. Phys. Chem. 34, 418 (1930). 
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CRYSTALLINE 


Parks? had indicated that such glucose glass 
preparations represent a mixture of the alpha- 
and beta-forms, and hence it seemed desirable to 
ensure here a final solution of mutarotated 
product comparable to that obtained when the 
crystals of a-glucose were previously dissolved 
in a mutarotating solvent. Our various results 
with the two glasses, as these were subjected to 
an aging treatment of from one to sixteen days, 
have been summarized in Table III. From these 
data it is noticeable that the heats of solution 
of Glass I averaged about 10 percent greater 
(numerically) than those of Glass IT. 

However, the aging process at 24°C had no 
appreciable effect on the character of our results, 
nor did the storage of samples of Glass II for 
thirteen days at 3° produce any significant change 
in its heat content as we measured this quantity. 
Previously, Parks and Reagh* had observed very 
appreciable changes in the viscosity of glucose 
threads when these were aged in the neighbor- 
hood of 20° to 30°; and more recently Nelson 
and Newton” have found that the average 
specific heat between 25° and 30° for a sample 
of glucose glass increased about 15 percent when 
stored at 25° for seven days. In view of our 
present results, it would seem that any appreci- 
able changes in the enthalpy of these glasses 
with aging must have developed largely within 
the first twenty-four hours. Moreover, our experi- 
mental method, which involved maintaining the 
glucose glass samples for one or two hours 
around 24° prior to a solution determination but 
following the storage for thirteen days at 3°, 
apparently precluded the observation of such a 
change of heat content as was noted by Nelson 
and Newton in their specific heat study. 


THERMODYNAMIC CALCULATIONS FOR THE 
CRYSTAL-TO-GLASS TRANSFORMATION 


The thermodynamics of the crystal-to-glass 
transformation have been previously considered 
by Parks, Snyder, and Cattoir,! but our new 
data now provide a more accurate basis for such 
a study with two somewhat different glass 
preparations. 

From the mean results for Glass I in Table III 
and for the crystals with mutarotation in Table II 


*E. W. Nelson and R. F. Newton, J. Am. Chem. Soc. 
63, 2178 (1941). 
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TABLE III. Heat of solution of glassy glucose samples 
in 0.0125N NH,OH at 25°C. 











Age of 

sample AH. for solution (cal. /g) 
Experiment No. in days* Glass I ylass IT 
GG-1 1 —6.15 
GG-2 1 —6.22 - 
GG-3 and 6 2 — 6.07 — 5.62 
GG-7 2 -—— — 5.62 
GG-8 2 — 58 
GG-4 3 —6.27 
GG-5 3 —6.21 
GG-9 11 —6.25 
GG-10 11 — 5.98 
GG-11 and 13 15 —6.00 —5.56* 
GG-12 and 14 15 — 6.08 — 5.62* 
GG-15 16 — — 5.61 
GG-16 16 —- —5.51 
Mean result —6.14 (+0.09) —5.57 (+0.05) 





* The aging of samples was carried out at a temperature of 24°C 
except in Experiments GG-13 and 14. These two samples of Glass II 
were aged for 13 days at 3° prior to the solution experiments. 


the enthalpy difference is 20.24 (+0.10) cal./g 
at 25°C. This yields for the molal process 


C gH 120¢ (a-crystals)>CeH 120¢ (glass) ; 
AH 9 = 3646 (+18) cal. (1) 


Above 25°C the glass softens and in the neighbor- 
hood of 40° the product in this process becomes 
obviously a very viscous liquid. An examination 
of the available specific heat data*” for crystal- 
line and undercooled liquid glucose indicates 
that the difference between the heat capacity of 
the undercooled liquid and the crystalline form 
is almost independent of the temperature; and 
accordingly we shall now take the value AC, = 35 
cal. per mole between 25° and 146° (the melting 
point of pure a-glucose crystals). Hence, by 
Kirchhoff’s law, we can calculate for the con- 
version of a mole of a-crystals into. liquid glucose 
at 146°C 


AHi419 = 3646+35 (146-25) 
= 7881 (+100) cal. (2) 


This result for the heat of fusion of a-glucose 
crystals is in fairly good agreement with the 
experimental one obtained by Thomas and 
Parks® with their radiation calorimeter. Their 
value for a glucose sample, somewhat contami- 
nated with decomposition products, was 41.7 cal. 
per gram at 141°, which corresponds to about 
7685 cal. per mole at 146°C. 

The entropy increase on fusion of the a-crystals 
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Taste IV. Thermodynamic data for the transformation 
a-glucose crystals—glassy glucose, per mole. 














Glass I Glass P-S-C 
AHos in calories 3646 3544 3893 
ASog in cal. /°K 6.87 6.63 7.46 
AF 93° in calories 1598 (+15) 1567(+415) 1669 (+50) 
AHj4i9 in calories 7881 7779 8128 
ASai9 in cal. /°K 18.80 18.56 19.39 
AF 4:9° in calories 0 


0 0 





at 146° follows directly from the data of Eq. (2): 
AS4i9 = 7881 /419 = 18.80 ©.%. (3) 


The entropy change forthe hypothetical con- 
version of a mole of the crystals into glucose 
glass at 25°C can be evaluated next with the 
aid of the AC, value used previously : 


= 6.87 (+0.07) e.u. (4) 


By means of the fundamental thermodynamic 
equation AF =AH—TAS, and the data of Eqs. 
(1) and (4), we can then obtain the molal free 
energy change for the crystal-to-glass transfor- 
mation at 25°C: 


AF o9s° = 1598 (+15) cal. (5) 


It may be noted that the relatively small un- 
certainty in this final free energy value is a 
consequence of a compensative effect in the 
corresponding uncertainties for the AH and 
ASog values. 

Similar data have also been worked out for 
the case of Glass II, and the preliminary results 
of Parks, Snyder, and Cattoir have now been 
revised with 35 for AC, instead of the value of 
33 calories per degree, which was originally used. 
The thermodynamic values for the three glasses 
are here summarized in Table IV to facilitate 
comparisons. 

From these tabulated data we obtain for the 
hypothetical transformation, per mole, 


Glass I1->Glass I], 
AHog = — 102 (+18) cal. ; 
ASoos = —().24 (+0.06) €.1; 
and 
AF a98° = —31 (+10) cal. 


Here the uncertainties are those which exist 


AND G. S. 


PARKS 


after the elimination of possible systematic 
errors, common to the data for both glasses and 
canceling in a transformation of this type. 

Such results strongly support the concept of 
a series of glassy states for a substance, with 
appreciably different enthalpy and free energy 
contents, as previously suggested by Parks and 
Reagh.* 


THE HEAT OF MUTAROTATION 
OF ALPHA-GLUCOSE 


Our data in Tables I and II provide new 
values for the enthalpy change when a mole of 
a-d-glucose, dissolved in water, is transformed 
into an equilibrium mixture of the alpha- and 
beta-forms. By use of the optical rotation values 
of Isbell and Pigman® we estimate that such an 
equilibrium solution contains 63.6 percent beta- 
and 36.4 percent alpha-glucose at 25°C. As our 
Pfanstiehl crystals initially appeared to be only 
about 96.0 (+1.0) percent a-glucose, we now 
calculate from Table I, for the mutarotation 
effect starting with a mole of pure a-glucose, 


ATI o93 = — (0.93) (180.156) (63.6/59.6) 
= —178 (+8) cal. 


Likewise, from the solution data for the pure 
alpha-crystals of Table II we obtain a somewhat 
more reliable value of 


AFT 3 = — (0.965) (180.156) = —174 (+5) cal. 


Accordingly, as a final, weighted mean we shall 
take AH/o3= —175 (+5) calories per mole. 

As noted in the previous discussion of the 
heats of solution of glucose crystals, several of 
the earlier investigators either overlooked this 
heat effect in mutarotation or assumed that it 
was zero. Brown and Pickering," however, ob- 
tained a value of 106 cal. per mole in 1897 and 
Sturtevant" reported 162 (+4) cal. in 1937. 
More recently, in a second paper® dealing: with 
precise calorimetry, he has found a result of 
174.2 (+0.9) cal., which agrees excellently with 
our final value of 175 (+5) calories for this 
mutarotation process. 


18H. T. Brown and S. Pickering, J. Chem. Soc. 71, 756 
(1897). 
4 J. M. Sturtevant, J. Am. Chem. Soc. 59, 1528 (1937). 
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An empirical relation is shown between contributions of various groups to the potential 
barriers of methyl groups and the minimum distance of approach of atoms of these groups 
to the methyl group. The use of this relation to derive barrier values for other molecules is 


demonstrated. 


I. INTRODUCTION 


HE existence of barriers hindering the rota- 

tion of groups about bonds has been known 

for some time. Methods have been described for 

the calculation of these potential barriers by the 

comparison of the results of thermodynamic 

measurements with the results of statistical 
mechanical calculations.' 

A method of estimating potential barriers 
where no thermodynamic data exist is highly 
desirable for the calculation of thermodynamic 
functions. The available evidence indicates that 
the equilibrium configuration of groups on op- 
posite ends of a CC axis is staggered rather than 
opposed.” A plausible explanation of the origin 
of these barriers is that they arise from repulsive 
forces between the atoms of the groups in rela- 
tive rotation. In the barriers considered here 
(methyl groups rotating on various molecular 
frames) these forces are presumably of the Van 
der Waals type, operating between the CH 
bonding electrons of the methyl group and the 
other atoms of the molecule.* 

It may be possible that energy vs. 
relationships, similar to the usual type of inter- 
atomic energy curves, exist for pairs of atoms on 


distance 


Harold Brunn Research Institute, 


* Present address: 
Mt. Zion Hospital, San Francisco, California. 
' K. S. Pitzer and W. D. Gwinn, J. Chem. Phys. 10, 428 


(1942). This is the most recent paper on the theory and 
gives references to earlier work. 

*'V. Schomaker and D. P. Stevenson, J. Chem. Phys. 8, 
637 (1940); K. W. F. Kohlrausch and H. Wittek, Zeits. f. 
physik,. Chemie B48, 177 (1941). 

* Effects due to unshared electron pairs have been 
ignored. The ultimate justification for this lies in the 
success achieved when only hydrogen-hydrogen interac- 
tions are considered. In other cases than those considered 
here it may be necessary to take into account the additional 
effect of dipole-dipole interactions: J. Y. Beach and D. P. 
Stevenson, J. Chem. Phys. 6, 635 (1938); J. Y. Beach and 
K. L. Palmer, J. Chem. Phys. 6, 639 (1938), 


opposite ends of a bond about which rotation is 
occurring. From such curves it would then be 
possible to calculate the energy changes which 
occur when groups rotate relative to one another. 
This approach was used in a recent paper by 
Aston, Isserow, Szasz, and Kennedy.‘ An inter- 
atomic potential of the form V=k/r;;" was as- 
sumed between pairs of atoms 77. The constants 
were evaluated from known barriérs in ethane 
and neopentane and other barriers were calcu- 
lated with fair success. A simpler and more 
empirical approach has been adopted by the 
present authors. This method reproduces the 
barriers of methyl groups in a certain class of 
molecules with somewhat greater accuracy and 
is easier to apply. 

As is the common practice, the assumption of 
independence of torsional coordinates is made, 
i.e., the absence in the potential energy function 
of cross terms involving products of torsional 
coordinates of different groups. Although such 
cross terms, particularly between adjacent tor- 
sional rotators, certainly do exist and have been 
evaluated for one simple case (see Pitzer, refer- 
ence 6), the difficulty of their evaluation and the 
fact that they would seem to be small compared 
to the “‘pure’’ potential terms have led to the use 
of effective barriers obtained by ignoring them. 


II. DEPENDENCE OF BARRIERS ON MINIMUM 
INTERNUCLEAR DISTANCES 


The barrier on a single methyl group is re- 
garded as being due to a summation of contribu- 
tions from interactions with groups on the other 
end of the rotational axis. For example, the 

1. & 


Kennedy, 


Aston, S. Isserow, G. J. Szasz, and R. M, 


J. Chem. Phys. 12, 336 (1944), 
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CONTRIBUTION TO V, colories 


Fic. 1. Relation of barrier contributions to minimum 
internuclear distance. 


barrier of 2750 calories found for ethane’ may 
be looked upon as three times the barrier offered 
to a methyl group by a single CH bond at the 
other end of the CC rotation axis. Hence, this 
barrier contribution may be taken as 920 cal. 
The observed barrier of 3300 cal. restricting the 
rotation of a methyl group in propane® may be 
regarded, similarly, as being a summation of two 
920-cal. contributions from the hydrogen atoms 
on the central carbon atom and a contribution 
due to the other methyl group. Thus this con- 
tribution has a value of 3300 —2 X 920 = 1460 cal. 
The observed barriers to the methyl group rota- 
tions in methanol and dimethyl ether may be 
regarded as being entirely due to the hydrogen 
atom and the methyl group, respectively, while 
those in neopentane and tetramethylsilane are 
three times the individual contribution from a 
methyl group on the other end of the rotation 
axis. 

The minimum distance of approach of a hy- 
drogen nucleus of the methyl group to any hy- 
drogen nucleus on the opposing group has been 
chosen as a convenient measure of the amount 
of interaction giving rise to a barrier contribu- 
tion. This minimum distance of approach is 
taken with respect to the rotation of the groups 
with all parameters other than rotational orienta- 
tion held constant. Thus for ethane the minimum 
distance is that from a hydrogen nucleus on one 


5G. B. Kistiakowsky, J. R. Lacher, and F. Stitt, J. 
Chem. Phys: 7, 289 (1939). 

6 G. B. Kistiakowsky and W. W. Rice, J. Chem. Phys. 8, 
610 (1940). A value of 3400 cal. was recently derived by 
K. S. Pitzer, J. Chem, Phys. 12, 310 (1944), 


FRENCH AND 


R. S. RASMUSSEN 


end to the corresponding hydrogen nucleus on 
the other end of the CC axis when the groups 
are in the opposed (D3,) configuration. For 
propane the minimum distance between methy| 
groups is calculated with both methyl groups in 
the “opposed” configuration with respect to the 
CHz group. Other molecules are treated in a 
similar fashion. When these minimum distances 
are plotted against the barrier contributions from 
the respective groups, the points may be fitted 
remarkably well by a single smooth curve (Fig. 
1). It is an oversimplification to assert that all 
the forces giving rise to the barrier are functions 
of this single parameter; nevertheless the excel- 
lence of the fit and the adequacy of the predic- 
tions made from the curve indicate that this 
parameter (minimum distance of approach) is a 
useful measure of the magnitude of the interac- 
tions of the groups involved. 

The barriers used in obtaining data for the 
construction of Fig. 1 are listed in Table I. All 
values given are taken directly from the indi- 
cated literature references with the exception of 
the one for dimethyl] ether, which together with 
methanol requires special discussion. 

For methanol an analysis of the fine structure 
of the infra-red absorption spectrum by Koehler 
and Dennison’ led to an estimate of 1340 calories 
for the barrier. Crawford® calculated a barrier 
of 3400 calories using the third law entropy ob- 
tained for liquid methanol by Kelley? and the 


TABLE I, Barrier contributions used in 
construction of Fig. 1. 











Barrier 
Measured contribu- 
barrier, tions, 

Compound cal. Distances, A cal. Reterences 
Ethane 2750 2.270 920 5 
Propane 3300 2.261, 1.991 930, 1440 6,4 
Neopentane 4800 1.923 1570 b 
Methanol 1350 2.030 1350 7 
Tetramethylsilane 1100-1500 2.556 370-500 c 
Dimethy! ether 700 1.774 2700 see text 








2 The values shown here for barrier contributions differ very slightly 
from those used in the text. This is due to a slight adjustment made 
in the methyl-to-methylene closest hydrogen distance, resulting from 
the non-tetrahedral angles about the central carbon, and making this 
distance somewhat shorter than the analogous ethane distance. 

+ J. G. Aston and G. H. Messerly, J. Am. Chem. Soc. 58, 2354 (1936) ; 
K. S. Pitzer, J. Chem. Phys. 5, 473 (1937); K. S. Pitzer, Chem. Rev. 
27, 39 (1940). 

¢J. G. Aston, R. M. Kennedy, and G. H. Messerly, J. Am. Chem. 
Soc. 63, 2343 (1941). 


7J. S. Koehler and D. M. Dennison, Phys. Rev. 57, 
1006 (1940). 

8B. L. Crawford, Jr., J. Chem. Phys. 8, 744 (1940). 

* K. K. Kelley, J. Am. Chem, Soc. 51, 181 (1929). 
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latent heat of vaporization measurements of 
Fiock, Ginnings, and Holton.!° However, since 
methanol has only a single hydrogen-methy] in- 
teraction at each maximum, it would be expected 
to have a barrier much smaller than that of 
ethane where three interactions contribute to 
each maximum. Much higher latent heats than 
those of Fiock, Ginnings, and Holton, which are 
a priort as acceptable, have been reported in the 
literature."! Third law entropies calculated using 
these data give a barrier in close agreement with 
the spectroscopic value of 1340 calories. For 
these calculations the frequency assignment of 
Noether™ and the moments of inertia shown in 
Table II were used. Consequently, the spectro- 
scopic value was adopted. Since there has been 
considerable lack of agreement on the barrier 
for this molecule, the data referred to above are 
discussed in detail in an appendix to this paper. 

For dimethyl ether an accurate third law en- 
tropy has been determined by Kennedy, Sagen- 
kahn, and Aston.'* Using the frequency assign- 
ment of Crawford and Joyce they obtained a 
barrier of 3100 cal. Unpublished infra-red data 
from these laboratories led to a frequency assign- 
ment sufficiently different from that of Crawford 
and Joyce to warrant recalculation of the barrier, 
leading to a value of 2700 cal. This frequency 
assignment together with the moments of inertia 
used is shown in Table II. The heat capacity 
data of Kistiakowsky and Rice on dimethyl 
ether vapor unfortunately cannot be used to 
determine the barrier, since contributions from 
the internal rotations in the temperature range 
of the measurements are nearly constant over a 
fairly broad range of barrier values. Small errors 
in either the frequency assignment or the heat 
capacity measurements would lead to large errors 
in the barrier. However, the calculated and ob- 
served heat capacities are in substantial agree- 





“ E. F. Fiock, D. C. Ginnings, and W. B. Holton, Bur. 
Stand. J. Research 6, 881 (1931). 

“E. Bartoszewicz, Rocznicki Chem. 11, 90 (1931); 
Chem. Abs. 25, 5343 (1931); S. Young, Proc. Roy. Dublin 
Soc. 12, 374 (1911). 

” H. D. Noether, J. Chem. Phys. 10, 693 (1942). 

*R. M. Kennedy, M. Sagenkahn, and J. G. Aston, 
J. Am. Chem. Soc. 63, 2267 (1941). 

*B. L. Crawford, Jr. and L. Joyce, J. Chem. Phys. 7, 
307 (1939). 

 G. B. Kistiakowsky and W. W. Rice, J. Chem. Phys. 
8, 618 (1940), , 


TaBLeE II. Moments of inertia and vibrational frequencies 
used in statistical calculations. 














Reduced 
Product of moment 
moments for internal Vibrational 
of inertia, rotation, frequencies, 
Compound g? cmé g cm? cm~! 
Dimethy! ether 1.708 -10-115 4.33 -10~40 412, 920, 1102, 
1180(2), 1260(2), 
1425(2), 1450(2), 
1462(2), 2930(6) 
Methanol 7.748 -107117 1.08 -10~40 1034, 1209, 1260, 
1340, 1430, 1455, 
1477, 2845, 2978(2), 
3683 
Methylamine 1.101 -10-"6 1.82 -10~4° 620, 783, 1064, 
1120, 1174, 1426, 
1460(2), 1625, 
2930(3), 3350(2) 
Acetone 1.398 -10714 4.99 -10~40 391, 488, 531, 787, 


900(2), 1066, 1070, 
1220, 1350, 1353, 
1420, 1423, 1430 
1438, 1710, 2930(6) 


7.867 -10 7115 4.90 -10~40 Same as those of Os- 
borne, Doescher, 
and Yost.¢ 


Dimethy! sulfide 











All physical constants used were taken from R. T. Birge, Rev. Mod. 
Phys. 13, 233 (1941). Thermodynamic contributions for restricted 
rotators were from the tables of Pitzer and Gwinn, reference 1. 

« Osborne, Doescher, and Yost, J. Am. Chem. Soc. 64, 169 (1942). 


ment. Consequently the value of 2700 calories 
was accepted. 

The molecular constants used in calculating 
the minimum internuclear distances used in 
Fig. 1 are given in Table III. The CH bond length 
of 1.093A from spectroscopic work on methane™® 
and the CC bond distance of 1.542A from x-ray 
diffraction work on diamond" were used through- 
out. The bond angles in methyl groups were as- 
sumed to be tetrahedral. All other values used 
and their sources are indicated in Table III. 
Interatomic distances were carried to 0.001A in 
all cases. This is well beyond the accuracy with 
which these parameters are known for most 
molecules. However, since small changes in these 
distances were significant, this procedure gave 
a greater degree of internal consistency than 
would have been obtained with values rounded 
off to the experimentally significant limit of 
0.01A. 


Ill. TESTS OF THE METHOD 


A number of calculations have been made to 
illustrate the use of the curve of Fig. 1 in esti- 


16 PD. M. Dennison, Rev. Mod. Phys. 12, 210 (1940). 

17 Value given by L. Pauling, Nature of the Chemical 
Bond (Cornell University Press, Ithaca, 1939), first edition, 
p. 150, 
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TaBLe III. Molecular structural parameters used in cal- 
culating minimum distances of approach of rotating 
groups. 











Molecule Parameters (lengths in A) 
Ethane standard 2 (see text) 
Propane <CCC, 111° 
Neopentane standard ies <CCC tetrahedral 
Methanol CO, 1.420; OH, 0.959; <COH, 104° 40’ 


CSi, 1.93; <CSiC tetrahedral 


Tetramethylsilane 
CO, 1.420; <COC, 111° 


Dimethyl ether 


Methylamine CN, 1.470; NH, 1.016; <CNH, 108°; 
<HNH, 108° 

Dimethylamine CN, 1.470; _ 1.016; <CNH, 108°; 
<CNC, 

Trimethylamine CN, 1.470; 'ECNC, 110° 

Methyl mercaptan CS, 1.823 SH, 1.35; <CSH, see text, 
Section III 

Dimethyl} sulfide CS, 1.82; <CSC, 96° 

Acetone CO, 1.24: CL, 1.323. <CCC, 120° 

Isobutane <CCC, 111° 30’ 

Trimethylborane BC, 1.56; <CBC, 120° 


PC, 1.87; <CPC, 100° — 

CSe, 1.93; SeH, 1.49; <CSeH, 90° 
CTe, 2.10; TeH, 1.68; <CTeH, 90° 
CGe, 1.98; <CGeC tetrahedral 
CSn, 2.18; <CSnC tetrahedral 
CPb, 2.23; <CPbC tetrahedral 


Trimethylphosphine 
Methy! hydroselenide 
Methy! hydrotelluride 
Tetramethylgermane 
Tetramethyltin 
Tetramethyllead 











These parameters were taken from the electron diffraction results 
summarized by: L. Pauling, Nature of the Chemical Bond (Cornell 
University Press, Ithaca, 1939), first edition, Chaps. III and V. V. Scho- 
maker and D. P. Stevenson, J. Am. Chem. Soc. 63, 37 (1941). 

In addition, spectroscopic results were utilized for certain of the 
parameters of methanol (OH distance; see A. Borden and E. F. Barker, 
J. Chem. Phys. 6, 553 (1939)), the amines (NH distance and HNH 
angle taken from NHs3), and methyl mercaptan (SH distance taken 
from HS). 


mating barriers and to demonstrate the accuracy 
and practicality of this method. The compari- 
sons of estimated barriers and those reported in 
the literature, for molecules other than those 
used in constructing the curve, are given in 
Table IV. In addition the potential barriers have 
been estimated for a number of molecules that 
are amenable to statistical treatment and for 
which no calorimetric measurements have been 
published. Molecular constants used in cal- 
culating minimum internuclear distances for 
these molecules are given in Table III. 

The following discussion of the check calcula- 
tions will illustrate the method as applied to the 
methylamines, methyl mercaptan, dimethy] sul- 
fide, acetone, and isobutane. 

The closest distance of approach of a methyl- 
group hydrogen to the hydrogen of an NH group 
linked to the methyl is 2.151A, using the mo- 
lecular dimensions in Table III for methylamine. 
Therefore, from Fig. 1, the barrier contribution 
of the NH bond is 1110 cal. The closest distance 
between a methyl-group hydrogen and a hy- 
drogen of another methyl-group bonded to the 
same nitrogen atom is 1.823A, leading to a 
contribution of 2130 cal. offered by the ~NCH; 
group to a methyl group in di- and: trimethyl- 
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amine. A small expansion of the CNC angle 
over the HNH angle was assumed for di- and 
trimethylamine analogous to the known in- 
crease in the CCC angle of hydrocarbons over 
the tetrahedral value. Summing the contribu- 
tions from the proper groups, the following bar- 
riers are obtained for methyl groups in the 
methylamines: 


monomethylamine = 2 X 1110 = 2220 calories, 
dimethylamine = 1130+ 2130, = 3260 calories, 
trimethylamine = 2 X 2130 = 4260 calories. 


These values are in close agreement with ex- 
perimentally determined values, as shown in 
Table 1V.18 

The third law entropy of methyl mercaptan 
leads to a barrier of 1460 calories. The value 
from Fig. 1 is 930-1010 depending on the CSH 
angle assumed. In the case of dimethyl sulfide, 
while an accurate third law entropy has been 
determined, the uncertainty in the CSC angle 
makes an accurate determination of the barrier 
impossible. Using an angle of 96° for calculation 
of the moments of inertia leads to a barrier of 
1990 calories for each methyl group. This same 
angle gives a minimum hydrogen-hydrogen dis- 
tance between methyl groups of 1.867A and a 
barrier of 1870 calories from Fig. 1, in satisfac- 
tory agreement. 

In the case of acetone, the third law entropy 
of the liquid,!® various reported heats of vaporiza- 
tion, the heat capacity of the gas,” and the 
frequency assignment shown in Table II gave 
barriers ranging from 560 to 1240 calories com- 
pared to 970 calories from Fig. 1 

For isobutane the third law entropy led to a 
barrier of 3870 calories.” Using a CCC angle of 
111°30' a _— of 3850 calories was obtained 
from Fi ig. 


‘Us sing the gaseous heat capacity data of Felsing and 
ime J. Am. Chem. Soc. 55, 4418 (1933), and the 
frequency assignment for methylamine given in Table II, 
values of C,(calc.)—C,(obs.) of +.59 at 0°C, —.05 at 
25°C, and —.52 at 50°C were obtained. A similar trend 

was noticed by Aston, Eidinoff, and Forster, J. Am. Chem. 
Soc. 61, 1539 (1939), in comparing calculated and observed 
Cosco Felsing and Jesson for dimethylamine. These trends 
cannot be accounted for by any reasonable modifications 
of barrier or frequency assignments. 

19K. K. Kelley, J. Am. Chem. Soc. 51, 1145 (1929). 

20K. Bennewitz and W. Rossner, Zeits. f. physik. 
Chemie B39, 126 (1938). 

217. G. Aston, R. M. Kennedy, and S, C, Schumann, 
J. Am. Chem. Soc. 62, 2059 (1940), 
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POTENTIAL BARRIERS 


IV. APPLICATION TO OTHER CLASSES 
OF MOLECULES 


This method for estimating barriers is, as 
emphasized above, applicable only to methyl 
groups, since it is only by virtue of the threefold 
symmetry that it becomes legitimate to split the 
barrier up into independent contributions from 
the groups giving rise to the barrier. However, 
some fairly successful attempts have been made 
to estimate barriers for other simple groups 
bearing hydrogens. 

Thus in ethanol the methyl-group barrier was 
calculated from the interactions with two methyl- 
ene-group hydrogens at 2.270A and one-third of 
the contribution from the hydroxyl hydrogen, 
whose closest approach is at 1.798A. Similarly, 
the barrier for the OH. group was calculated from 
Fig. 1, using an average of values obtained from 
the minimum distances of approach of the OH 
hydrogen to the methylene-group hydrogens and 
to the nearest methyl-group hydrogen : 

2400+ 2 X 1340 


Von = = 1690 cal. 
3 





2400 
Vue = 2 X916+— = 2630 cal. 
3 


In isopropyl alcohol the analogous values are 


TABLE IV. Barriers derived from Fig. 1, and comparison with experimental values. 
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Vou =2070, Vwe=3240, and in tertiary butyl 
alcohol Von = 3400 and Vy. = 4010. These values 
for the OH group barriers are believed to be 
reliable. If, however, the oxygen atom has any 
large effect the methyl barrier estimates may be 
low. Preliminary comparisons of statistical cal- 
culations for these alcohols with the results of 
all known measurements are satisfactory. The 
data on these alcohols are from various sources, 
which are often in serious disagreement. Hence, 
a complete treatment of the data, which will be 
presented in a later paper, is required to clarify 
the picture. 

In molecules such as butane where the rotating 
groups are highly unsymmetrical the present 
method does not seem to be applicable, and in 
such cases it will probably be necessary to use a 
more difficult but more fundamental treatment 
such as that of Aston. 

A consideration of the barriers hindering the 
rotation of methyl groups attached to ethylenic 
frames has been omitted from this correlation 
for several reasons. In the rotation of these groups 
there may be six maxima per rotation and the 
barrier contributions might be expected to be 
different from those predicted from Fig. 1, which 
was constructed from data on molecules having 
three maxima per rotation. 











Total barrier 








Contributions from Measured 
from Fig. contributions, barrier, 

Compound Distances, A cal. cal. cal. References 
Methylamine 2.151 1110 2220 1520 or more a 
Dimethylamine 2.140, 1.823 1130, 2130 3260 3460 b 
Trimethylamine 1.823 2130 4260 4270 c 
Methyl mercaptan 2.208-2.262* 1010-930 1010-930 1460 d 
Dimethyl sulfide 1.867** 1870 1870 2000 e 
Acetone 2.271 920 920*** 560-1240 see text, 

Section ITI 
Isobutane 2.232, 1.991 970, 1440 3850 3870 f 
Trimethylborane 2.302 860 1720 
lrimethylphosphine 2.098 1210 2420 
Methyl hydroselenide 2.340 800 800 
Methyl hydrotelluride 2.548 470 470 
Tetramethylgermane 2.638 ~400 ~1200 
letramethyltin 2.965 ~0 ~~ 
Tetramethyllead 3.046 - ~0 ~0 














* Assuming a CSH angle of 90° and 92° 20’, respectively. 
* Assuming a CSC angle of 96°. 


*** Assuming the whole barrier to be due to the other methyl group. The barrier derived on this assumption fits the available experimental data 


within their probable errors. 


- J. G. Aston, C. W. Siller, and G. H. Messerly, J. Am. Chem. Soc. 59, 1743 (1937); J. G. Aston and P. M. Doty, J. Chem. Phys. 8, 743 (1940). 
J. G. Aston, M. L. Eidinoff, and W. S. Forster, J. Am. Chem. Soc. 61, 1539 (1939). 

‘J: G. Aston, M. L. Sagenkahn, G. J. Szasz, G. W. Moessen, and H. F. Zuhr, J. Am. Chem. Soc. 66, 1171 (1944). 
H. Russell, Jr., D. W. Osborne, and D. M. Yost, J. Am. Chem. Soc. 64, 165 (1942). 

*D. W. Osborne, R. N. Doescher, and D. M. Yost, J. Am. Chem. Soc. 64, 169 (1942). 
J. G. Aston, R. M. Kennedy, and S. C, Schumann, J. Am. Chem. Soc. 62, 2059 (1940). 
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V. CONCLUDING REMARKS 


This method, like any other method of esti- 
mating barriers based on an exact knowledge of 
interatomic distances, is very sensitive to changes 
in these parameters, particularly in the region of 
short internuclear distances. For this reason an 
effort has been made to use a consistent set of 
molecular constants. Values consistent with those 
which entered into the construction of Fig. 1 
must be used in making barrier estimates from 
this curve. The arbitrariness in the choice of 
molecular constants is particularly evident in 
the choice of hydrogen angles. It is impossible, 
at present, to state the precise values for these 
angles in complex molecules because they cannot 
be obtained by current spectroscopic and elec- 
tron diffraction methods. It has been assumed, 
for example, that the COH angle is the same as 
the water angle and that HCC and HCH angles 
are tetrahedral, except in cases where this is 
obviously untrue (such as the HCC angle of the 
tertiary hydrogen in isobutane). 

An additional point of question lies in the 
uncertainties of the values of the barriers used 
in constructing Fig. 1. These uncertainties may, 
in some cases, be as high as +500 calories. How- 
ever, the excellence of fit of the barrier contribu- 
tions to a smooth curve, and the relatively small 
deviations between the barriers estimated from 
the curve and those estimated by the usual 
methods from calorimetric data indicate that 
this method will yield values which are entirely 
adequate for practical thermodynamic calcula- 
tions. 

The concept that repulsive forces give rise to 
potential barriers in the types of molecules con- 
sidered is supported further by the electron dif- 
fraction evidence that CCC angles are greater 
than tetrahedral in propane, isobutane, and ter- 
tiary butyl chloride.” In addition the change in 
some methyl rocking frequencies from approxi- 
mately 1160 wave numbers in methyl groups at- 
tached to secondary and tertiary C atoms to 
aL, Pauling and L. O. Brockway, J. Am. Chem. Soc. 
59, 1223 (1937); J. Y. Beach and D. P. Stevenson, J. Am. 


Chem. Soc. 60, 475 (1938); J. Y. Beach and J. Walter, 
J. Chem. Phys. 8, 303 (1940). 
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approximately 1250 in methyl groups attached to 
quaternary C atoms is readily interpreted as 
being due to the increased repulsion from ad- 
jacent groups. 


APPENDIX: THE BARRIER IN METHANOL 


The entropy of methanol vapor at 298.16°K 
and one atmosphere can be derived from avail- 
able data as follows: 


Entropy of liquid, 298.16°K 
(Kelley?) 

Entropy of vaporization! !! 

Gas imperfection correction 
(Berthelot Eq.) 0.02 

Entropy change, vapor pressure 
to one atmos. 


30.26 e.u. 
30.015, 30.447, 30.688 


— 3.586 


Entropy of vapor at 298.16°K 


and one atmos. 56.71, 57.14, 57.38 e.u. 


The three values given for entropy of vaporiza- 
tion are those derived from the latent heat 
values of Fiock, Ginnings, and Holton,!° Young," 
and Bartoszewicz,'! respectively. These values 
are 8949.4 cal., 9078 cal., and 9150 cal. 

The statistically calculated entropy, except for 
the internal rotational contribution, is (using 
constants and frequencies shown in Table II): 





Vibrational 0.220 e.u. 
Translational-rotational 55.291 
55.511 e.u. 


Thus the contribution from internal rotation is 
1.20, 1.63, or 1.87 e.u., respectively, obtained 
from the three latent heat values. The first value, 
1.20 e.u., derived from the Fiock, Ginnings, and 
Holton datum, leads to a barrier of 2360 cal. 
The spectroscopic value for the barrier (1340 
cal.) gives an internal rotational contribution of 
1.65 e.u., in essential agreement with the latter 
two values. The discrepancy between the barrier 
(2360 cal.) derived here using the Fiock, Gin- 
nings, and Holton latent heat and that given by 
Crawford® (3400 cal.) derived from the sane 
latent heat and liquid entropy is due mainly to 
the more recent physical constants and frequency 
values used here. 
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The entropies of the formic and acetic acid dimers are used to find a minimum value for 
the force constant of the association bond. With the aid of a normal coordinate treatment the 
probable value is found to be 4X 104 dynes/cm. This agrees with a rough estimate from the heat 
of dimerization. A Raman band in the formic acid spectrum at 200 cm~! may be due to the 


association vibrations. 





N earlier papers,! the calculation of reasonably 

accurate entropies from vapor densities, vapor 
pressures, and thermal data has been described 
for the momomers and dimers of formic acid and 
acetic acid. The dimer entropies, for the hypo- 
thetical perfect gases at 25°C and atmospheric 
pressure, were found to be 83.1+0.1 and 101.0 
+0.3 e.u., respectively, for the two acids. 

These entropies may be used to estimate the 
stretching force constant of the hydrogen bonds 
by which the two monomer units are thought to 
be held together in the dimer. The procedure is 
to find the entropy and frequency of the sym- 
metrical vibration of extension and compression 
of the association bonds (»; in Fig. 1) and from 
this frequency to calculate the force constant by 
means of the familiar relation 


2K =49'*vy=27'M, (1) 


where M is the mass of the monomer. 

From the total dimer entropy the contribu- 
tions of translation and rotation are subtracted 
to give the internal entropy. It is then assumed 
that the internal entropy of the monomer is not 
appreciably changed by combination into the 
dimer, so that the dimer internal entropy is the 
sum of twice the monomer internal entropy plus 
the entropy of six vibrations dependent upon the 
strength of the association bonds. These vibra- 
tions are illustrated in Fig. 1. For all the mono- 
mer frequencies except the rotational motion of 
the hydroxyl group this assumption is accept- 
able without comment. Extension to include the 
torsional motion follows at once from the already 
low entropy (0.75 in formic acid) which could 
only be further lowered to a minor extent by the 
restrictive forces resulting from the dimerization. 


(1983) Halford, J. Chem, Phys. 9, 859 (1941); 10, 582 


Table I shows the quantities used to obtain 
the entropy of the association frequencies. The 
monomer internal entropy of both acids appears 
in the first two rows, with the value for hydroxyl 
torsion set in each case at 0.75, as originally cal- 
culated for formic acid, because of the relatively 
high probable error in the corresponding quantity 
(1.2) in the acetic acid calculation. This choice 
has no effect on the conclusions reached in this 
paper. The first four rows in the dimer summa- 
tion in the lower part of the table contain the 
translational and external rotational entropy, 
with the constant adjusted to permit direct sub- 
stitution of the monomer weight in grams for 
M and of 10** times the principal moments of 
inertia in c.g.s. units for A, B, and C. @ is the 
symmetry number. 

A model based upon the original structure 
proposed by Pauling and Brockway,’ rather than 
the newer and less symmetrical one found by 














Fic. 1. Vibrations of the formic acid dimer. 


*L. Pauling and L. O. Brockway, Proc. Nat. Acad. Sci, 
20, 336 (1934). 
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Fic. 2. Coordinates and distances in the formic acid dimer. 


Karle and Brockway,’ is used to calculate the 
principal moments of inertia. The rotational 
entropy does not vary appreciably with minor 
changes in the model. Dimensions used are 
(H—C)=1.09, (C—O)=(C=O)=1.28, and 
(C—C) =1.54A units, with the O—C—O angle 
set at the rounded value of 120°. The oxygen 
atoms of the individual monomer units are taken 
to be linked through the bonding hydrogen at a 
distance of 2.67A units, with the hydrogen atom 
placed at 1.00 unit from the nearer oxygen. For 
the principal axes it is accurate enough to take 
the common bisector of the two O—C—O angles 
and its perpendiculars in and at right angles to 
the molecular plane. The two axes in the molec- 
ular plane are thus placed slightly out of their 
true positions because of the location of the 
hydroxyl hydrogens. 

By analogy to the properties of ordinary 
valence bonds, »; and v3 should represent con- 
siderably higher frequencies than the remaining 
four vibrations. It is safe to assume therefore 
that the entropy of the first frequency is less 
than one-sixth of the total shown in Table I for 
the association vibrations. Thus S; is less than 
2.53 e.u. for formic acid, and, from the Einstein 
function, v; is greater than 162 cm. The corre- 
sponding value for K is 1.7X10* dynes/cm. 
Similarly, for acetic acid, S; is less than 3.09, 
is greater than 121, and K is above 1.1X10*. 
The minimum is therefore determined by the 
more accurate formic acid entropy. 

The probable value is considerably larger be- 
cause at least one of the bending motions would 
be expected to have a much larger entropy, 
leaving less for »;. It is more difficult to estimate 
a maximum. From symmetry considerations 7 


8 J. Karle and L. O. Brockway, J. Am. Chem. Soc. 66, 
574 (1944). 
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should be active in Raman scattering but there 


is no way of predicting how intense the effect 
should be. It seems improbable, even if the 
intensity is low, that the Raman line would 
have escaped detection if the frequency were 
greater than 250 cm~, for which K for formic 
acid is 4.210‘. For acetic acid the constant 
derived from this frequency would be higher. 

The result of this estimate, K =3.0+1.3 X 10', 
can be supported in two independent ways, and, 
in addition, by a more detailed interpretation of 
the entropy by means of the normal coordinate 
method. A Raman band has been reported for 
liquid formic acid at 180+135 by Edsall,‘ at 200 
by Kohlrausch, Koppl, and Pongratz,®’ and at 
197 cm by Morino and Mizushima.® According 
to Hibben’ this may be analogous to a similar 
scattering observed for water, for which he sug- 
gests as a possible cause a hindered rotation in 
the liquid state. It is also possible that the effect 
may be due to the vibrations under discussion 
here. If so, a frequency of 200 cm would indi- 
cate for the hydrogen bond a stretching force 
constant of 2.7 10' dynes/cm. 

A second independent supporting estimate is 
based upon the assumption that the ratio of 
force constant to energy of dissociation is about 
the same for the association bond as for ordinary 
hydrogen valence bonds. For the CH, OH, CIH, 


TABLE I. Entropy of association vibrations at 25°C. 











Formic acid Acetic acid 
Monomer 
Vibration 0.60 3.60 
Hydroxyl torsion PY 5 Ri 
Methy] torsion a= 2.99 
Internal entropy 1.35 7.34 
Dimer 
Constant 49.87 49.87 
3Rin2M 13.48 14.27 
$R In ABC 3.22 5.04 
—Rine — 1.38 — 1.38 
Internal, monomer 2.70 14.68 
Association bonds 15.21 18.56 


83.10 101.04 








‘i T. Edsall, J. Chem. Phys. 4, 1 (1937). 
. W. F. Kohlrausch, F, Koppl, and A. Pongratz, 
zeits < physik. Chemie 21B, 242 (1933). 
6 Y. Morino and S. Mizushima, Sci. Pap. Inst. Phys. 
Chem. oe (Tokyo) 32, 33 (1937). 
7J. H. Hibben, The Raman Effect and Its Chemical 
A pplications (Reinhold Publishing Corporation, New York, 
1939), p. 180. 
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TABLE II. Molecular constants of monomer units. 


TABLE III. Equations for association frequencies. 








Formic acid Acetic acid 


10-22X (formic acid) 10%) (acetic acid) 





1072M 0.7593 0.9906 
10*D 0.7906 1.6053 
10**E 0.6685 0.7214 
10°*F 0.1221 0.9365 
10a 1.109 1.109 
10% 1.183 1.675 
10% 0.183 0.675 
108d 1.518 2.010 








and BrH bonds this ratio has an average value 
near 5.7 dynes/cm/cal. The heat of dissociation 
of the formic acid dimer, 7150 cal. per hydrogen 
bond, is taken to be near enough to the energy of 
dissociation and is multiplied by the above ratio 
to give K=4.1X10*. 

Finally, to provide support for the assumption 
that the entropy of v; should be considerably less 
than the average for the six association fre- 
quencies, a normal coordinate treatment is un- 
dertaken to determine what spread of frequencies 
is consistent with a simplified but reasonable 
potential function. In addition to its application 
to the present problem, this study of a planar 
system with a center of symmetry but no other 
symmetry elements is of interest for itself. 

The molecular model is shown in Fig. 2 which 
gives part of the symbolism including the co- 
ordinates applied to one monomer unit in setting 
up the initial equations. Because it has been as- 
sumed that the association vibrations are inde- 
pendent of other molecular motions, the model 
consists of two rigid bodies vibrating against 
each other. For such a system the kinetic energy 
may be conveniently expressed in terms of the 
translation and rotation of the individual rigid 
bodies. X, Y, and Z represent linear displace- 
ments of the centers of mass, while a, 8, and y 
are angular displacements about the principal axes. 

The kinetic energy is given by the following 
equation : 
2T = M(X2+ YP+Ze+X2+ Y2+Z,?) 

+D(ay?+ a2?) +E(BP+8)+F(y2 +727). (2) 


Subscripts distinguish the two monomer units. 

D, E, and F are moments of inertia about the 

axes PZ, PX, and PY, respectively, of Fig. 2. 
For the potential energy, 


2V = K[(x1—x2)?+(x3—x4)?] 
+2k(y2+y2+ye+y2 
+2P+22 +29 +20). (3) 


5.268K 4.038K 
2.497k 2.900k 
6.223K+3.625k 3.065K +4.063k 
21.941k 5.645k 
26.605k 9.593k 
4.228k 4.193k 








Here the subscripts indicate the terminal points 
of the association bonds, with even numbers 
applying to one monomer unit and odd numbers 
to the other, as shown in Fig. 2. To obtain Eq. 
(3) the general quadratic potential function has 
been simplified by eliminating all cross products 
except the x:x2 and x3x, terms and reducing the 
number of constants to two. For this form, K is 
evidently the stretching valence force constant, 
while k, which is concerned with lateral displace- 
ments, must be a bending force constant. Bend- 
ing constants defined in this way are about one- 
fifth to one-tenth as large as stretching constants, 
and may be converted to angular valence force 
constants by multiplying by the squares of the 
appropriate bond lengths. When, as in the pres- 
ent problem, the same constant is used for all 
the angular displacements, the factor required 
for this conversion is not evident, and the rela- 
tion to any specific angular valence force con- 
stant is not definite. The product k(b?+c) in 
Eqs. (7) and (8) should have about the right 
magnitude for an angular valence force constant 
because 6 and ¢ are of the correct order for inter- 
atomic distances. Equations (8) show the correct 
form for the angular vibration between like rota- 
tors with the reduced moment of inertia D/2, etc. 
These solutions were put into this form by intro- 
ducing the factor 2k rather than k into Eq. (3). 
Six of the twelve coordinates in the kinetic 
energy equation are eliminated by means of the 
following linear and angular momentum con- 
ditions: 
¥ =X, = — - 
Y=Y,=-Y, 2MdY =D(a,+d2) (4) 
The substitutions a;=6,;+62, a2=6:—62, yi=¢1 
+. and y2=¢1—¢2 are then made, leading to 
the expression 
2T =2MX?+2D(1+D/Ma?)62+2D6.2 
+2EB+2F(1+F/Md)¢r+2F¢?. (5) 


B=61= — Be 





‘: 


TABLE IV. Trial solutions. 
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Formic acid, K =5k 


v/vi vy (cm~) S298 
1.000 229 1.89 
0.308 71 4.13 
1.148 263 1.65 
0.913 209 2.06 
1.005 229 1.89 
0.401 92 3.62 

K=3.5X 104 15.24 

Formic acid, K = 10k 

1.000 293 1.46 
0.218 64 4.35 
1.118 328 1.28 
0.645 189 2.24 
0.711 208 2.06 
0.283 83 3.82 

K =5.7X 10" 15.21 

Acetic acid, K =5k 

1.000 193 23% 
0.379 73 4.07 
0.980 189 2.24 
0.529 102 3.42 
0.689 133 2.91 
0.456 88 3.70 

K =3.3X 10 18.55 

Acetic acid, K = 10k 

1.000 246 1.76 
0.268 66 4.27 
0.927 228 1.89 
0.374 92 3.61 
0.482 120 3.11 
0.322 79 3.92 

K=5.3X 104 18.56 








For the transformation of the potential energy 


the required equations are 


x1 =X —day yi =ba, z,=ap+by1 
X2x=—X—da. Ye=—Cay - 22= —aB—cy2 (6) 
x3=X+aa, Y3=Cay 23= —aB+cy 

x4= —X+daz2 Ys= —CQ2 24=aB—byo, 


together with the definition equations for 6), 62, 
¢1, and ¢». After the transformation the potential 
energy has the form 
2V =8KX?+4k(? +c?) (62+ o°+ ¢2") 
+4[2Ka?+k(b?+c?) ]0.? 
+8ka(b—c)Bgy2+8ka’B’. (7) 
The solutions for \=47’v? for the first four 
vibrations of Fig. 1 take the forms 
\.=4K/M Ao = 2k(b?+c*?) /L[D'1+D/Ma?) ] 
As=[4Ka?+2k(b?+<c*) ]/D (8) 
Ag = 2k(b?+0?)/LF(1+F/Ma?) }. 















































For Xs and Xz, 


EX —4ka? 
‘ —2ka(b—c) 


r-ne+re)|~- © 
— 2k(b’+-c’) 

The masses, moments of inertia, and distances 
required for numerical solution of the above equa- 
tions are collected in Table II. The moments 
D, E, and F have been calculated with the aid 
of the dimensions listed earlier for the dimer and 
are subject to a similar but relatively larger error 
because of the deviation of the selected rotational 
axes from the principal axes. The error is again 
taken to be negligible. The distances a, b, and c 
are coordinates referred to the axes PX and PY, 
as shown in Fig. 2, and d is the distance between 
O, the center of mass of the dimer, and P, that 
of the monomer unit. 

Substitution of the constants of Table II into 
Eqs. (8) and (9) yields the simplified expressions 
of Table III. 

In principle the two potential constants could 
be found by trial from the two entropy totals of 
Table I, but the result so obtained would have 
little meaning because of the relatively large un- 
certainty in the acetic acid entropy. It seems 
better to assign a reasonable value to K/k, to 
evaluate from this the ratios \/A; and v/y, and 
to find by trial from the entropy-frequency curve 
a set of frequencies in the corresponding ratios 
which will yield the correct entropy. Table IV 
shows in successive columns the frequency ra- 
tios, the frequencies, and the entropy summation 
for four such trial solutions. 

The assumed ratios of 5 and 10 for K/k are 
of the order usually found for ordinary valence 
bonds. For each ratio the spread of the calculated 
frequencies is several-fold and »; appears at or 
near the maximum, accounting for considerably 
less than one-sixth of the entropy of the asso- 
ciation frequencies. It is therefore probable that 
v; for formic acid lies in the region between 200 
and 250 cm~ and that K is not far from 4X 10° 
dynes/cm. 

In the absence of a more accurate method of 
evaluation this provides a satisfactory approxi- 
mation to the stretching potential constant of the 
association bond in the carboxylic acid dimer. 
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New Units for the Measurement of Radioactivity 


E. U. ConDON AND L. F. Curtiss 


National Bureau of Standards, U. S. Department of Commerce, 
Washington, D.C. 


May 20, 1946 


T has become the custom to express the strength of 

radioactive sources in terms of curies. This is an 
erroneous use of this unit since by original definition the 
curie is that “amount of radon in equilibrium with one 
gram of radium”! as defined by the Radiology Congress 
in Brussels in 1910. Therefore, the curie can be used only 
to represent a rate of disintegration in the radium family. 
It then represents the disintegration rate of radium or its 
products in equilibrium. Such a use has been indorsed by 
the International Radium Commission. 

The quantity to be specified in designating the strength 
of radioactive sources in general is the disintegration rate, 
determined by the decay constant and the number of 
atoms of the radioactive isotope in the source. This is 
simply a number, and therefore to establish a suitable unit 
the only requirement is to select a convenient number of 
disintegrations per second and give it a name. In selecting 
this number, consideration should be given to insure that 
it can be readily expressed in sub-multiples and multiples 
by the usual prefixes, kilo, milli, micro, etc. A number 
which fits this requirement is 10°. Since the curie was 
named in honor of M. and Mme. Curie, the co-discoverers 
of radium, it is natural to select the name “rutherford” 
for the new unit. The appropriate abbreviation is “rd,’’ 
which conflicts with the abbreviation of no other well 
accepted physical unit. The micro-rutherford would be- 
come one disintegration per second, a convenient number 
to remember. Furthermore, the rutherford itself is a small 
unit of the order of magnitude of many sources used in 
laboratory measurements. It is sufficiently different in 
size from the curie that no confusion can arise with the 
curie in connection with measurements of activities in the 
radium family. Large sources would require the use of 
positive powers of ten, which would be preferable to the 
use of a large unit requiring negative powers of ten. 

It should be pointed out that the continued use of the 
curie for all radio-isotopes not only requires a redefinition 
of the curie, but in addition the value of the curie is un- 
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certain to at least 4 percent and values are in current use 
well outside this limit. The rutherford provides a definite 
unit. In addition to eliminating the undesirable use of the 
curie, the proposed unit also eliminates the basic necessity 
for measuring radio-isotopes in terms of a standard. Any 
measuring device which will determine the total number 
of disintegrations per second will provide directly the 
strength of the source in rutherfords. A counting arrange- 
ment for which the solid angle factor is known is an ex- 
ample. Radioactive standards may be used to determine 
this factor for a given geometrical arrangement, but other 
methods are also available. 

In the measurement of sources of gamma rays, the 
roentgen has gained increasing use, largely because this 
unit is independent of the quality (electron volts) of the 
gamma radiation. There is need for a unit in which the 
intensity of gamma-ray sources can be expressed to elimi- 
nate the use of the curie for this purpose. An obvious unit 
derived from the definition of the roentgen is a roentgen- 
per-hour at one meter. The roentgen-per-hour at one 
meter can be abbreviated r.h.m., which again is not 
readily confused with any other common abbreviation. 
It has been suggested that this abbreviation can be pro- 
nounced “rum.” It should be noted that a gamma-ray 
source equal to one r.h.m. will have a gamma-ray strength 
1.18 times that of 1 curie of radium. Therefore, the 
roentgen-per-hour at one meter has the same order of 
magnitude as the curie in the measurement of gamma-ray 
sources. 

The National Bureau of Standards, at the suggestion of 
the Committee on Radioactivity of the National Research 
Council, recommends the general use of these units. 

1E. Rutherford, ‘‘Radioactive Substances and Their Radioactivity” 


(1913), p. 479. ‘The Radioactive Constants as of 1930,’’ International 
Radium-Standards Commission Report, Rev. Mod. Phys. 3, 427 (1931). 





The Electron Affinity of Amide Radical and the 
Chemical Stability of Metal-Ammonia Solutions 


RICHARD A. OGG, JR. 
Department of Chemistry, Stanford University, California 
May 13, 1946 


EW experimental studies have been carried out, deal- 

ing with concentration and temperature dependence 
of electrical conductivity and of the near ultra-violet 
absorption spectrum of liquid ammonia solutions of sodium 
and potassium amides. These substantiate older observa- 
tions! that the amides behave as weak electrolytes, and 
indicate that the dissociation constants increase with rising 
temperature. They also show that the constituent respon- 
sible for the yellow color of such solutions is the amide ion, 
NH:;-. Of particular importance is the observation that 
light absorption by the amide ion results in a striking 
photo-conductivity. This fact supports the interpretation 
of the observed continuous absorption spectrum as an 
electron affinity spectrum,? corresponding to the physical 
process 

NH2-(NH3;)+hv—NH2(NH;)+e(NHs). 


No over-all photochemical reaction is observed, suggesting 
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that recapture of conduction electrons by the free amide 
radicals must be rapid. 

The long wave-length threshold of the amide ion light 
absorption is at about 4200A, corresponding to an energy 
of some 68 kilocalories per mole. This value would appear 
to represent a numerical upper limit for the energy change 
associated with the above dissociation process. Some re- 
lated considerations employing Hess’ law are of paramount 
importance in understanding the remarkable chemical 

stability of metal-ammonia solutions. In the following 
equations, the inequality signs refer to algebraic values of 
the corresponding reaction energies, given in kilocalories 
per mole: 


NH:;(g)—>NH)2(g) +H (g) E,<—84, (1) 
NH2(g)—>NH2(NHs) Ex~ +6, (2) 
NH; (1) ~NH:(g) E;~ —6, (3) 

H (g)—H(NHs) Ew 0, (4) 


e(NH3) +NH2(NH3)—NH2-(NHs) Es<+68, (5) 





Reaction (6), being endothermic by at least 16 kilocalories, 
must have an activation energy at least as great as this 
value. Since the dissociation energy of the “first” N—H 
bond in NH; is without doubt greater numerically than 
the average value* used for reaction (1), a more probable 
estimate of the minimum activation energy of (6) would 
be of the order of 30-40 kilocalories per mole. That such 
large values actually obtain is indicated by the author's 
observations that highly purified metal-ammonia solutions 
preserved in sealed vessels at room temperature require 
several months for practically complete discharge of the 
characteristic blue color. The probable role of catalysts is 
to lower the energy deficit of reaction (6) by adsorption of 
hydrogen atom. In support of this postulate is the fact 
that the metals most suitable for hydrogenation catalysts 
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are also the most effective in promoting amide formation 
in metal-ammonia solutions. 

If the far ultraviolet absorption spectrum of aqueous 
hydroxyl ion‘ is interpreted as an electron affinity spec- 
trum, entirely analogous considerations indicate the reac- 
tion 

e(aq) +H20(/)>+OH~ (ag) +H (aq) 
to be considerably exothermic and hence probably asso- 
ciated with a negligible activation energy. This affords a 
ready explanation of the impossibility of preparing aqueous 
metal solutions containing ‘‘solvated” electrons, and sug- 
gests that affinity spectra of aqueous negative ions should 
not be attended by photo-conductivity. 

If the solvation energy of conduction electrons in am- 
monia is estimated as less than 10 kilocalories,5 then the 
numerical upper limit for the sum of the electron affinity 
of amide radical and the solvation energy of amide ion in 
ammonia is roughly 68+6+10=84 kilocalories. Any rea- 
sonable estimate of the amide ion solvation energy results 
in an extremely small electron affinity for the amide 
radical—of the order of 20-30 kilocalories. Such a value 
would fall in logical sequence with the respective figures 
of 95 and 48 kilocalories for the electron affinities* of F 
and OH. So small a value for the electron affinity of amide 
radical, as compared with the ionization potentials of 
alkali metal atoms, would suggest that the binding in the 
alkali metal amides is predominantly covalent. The experi- 
mentally observed weak electrolyte behavior in liquid 
ammonia solution supports this view, as do the relatively 
high volatibility and low melting points of the pure sub- 
stances.® 


1C, A. Kraus, The Properties of Electrically Conducting Systems 
(Chemical Catalogue Company, 1922). 

2 J. Franck and G. Scheibe, Zeits. f. physik. Chemie A139, 22 (1928). 

30. K. Rice, Electronic Structure and Chemical Binding (McGraw- 
Hill Book Company, Inc., New York, 1940). 

4H. Ley and B. Arends, Zeits. f. physik. Chemie B6, 240 (1929). 

5R. A. Ogg, Jr., J. Chem. Phys. 14, 295 (1946). 

6 F, W. Bergstrom and W. C. Fernelius, Chem. Rev. 12, 43 (1933). 
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